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Abstract 


Nano-crystalline metals (NCs) exhibit radiation-tolerance due to the sink of grain boundaries (GBs) for 
radiation-induced defects such as self-interstitial atoms (SIAs) and vacancies (Vs). However, the relevant 
mechanisms for the radiation damage accumulation and GB structural relaxation under high radiation field 
in NCs are still not well understood due to the lack of self-consistent across-scale techniques for simulating 
radiation-induced microstructures evolution. In this article, by combining coarse-grained and atomistic 
simulations, we proposed a coupling method to investigate the evolution of the microstructure and 
SIA/V-GB interaction under cumulative irradiation in NC iron. The SIA overloaded effect was revealed in 
iron GBs at a high radiation dose rate and/or low temperature. Two types of GB structural response were 
observed to cumulative irradiation. With the SIA accumulated at the GB, the new GB phase formed and then 
a critical concentration of the SIA at the GB transited to the small quantity of the V during the GB structural 
recovery, accompanied by the local GB motion. Consequently, the GB's role for Vs nearby alternated 


between the trapping and annihilation center with radiation dose. Alternatively, GB developed to a 


disordered structure after trapping abundant SIAs. The GB response pattern to cumulative irradiation that is 
related to the SIA formation energy at the GB or the GB thermal stability is well manifested in the 
cumulative distribution function of the defects formation energy and its energy level density. The present 


work reveals the dynamic healing picture for radiation damage near the GB under cumulative irradiation. 
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1. Introduction 

Current and future generation nuclear reactors call for structural materials with outstanding properties 
that can withstand hundreds of atomic displacements per atom (dpa) produced by neutron damage []. 
Radiation to the material creates defects of self-interstitial atoms (SIAs), vacancies (Vs) and their clusters 
(V,/SIA,) in the material []. These defects often induce the degradation of materials performance in terms of 
irradiation dislocation hardening, creep, void swelling, and helium embrittlement []. Engineering materials 
with a high density of defect sinks, such as GBs in nano-crystalline metals (NCs) with the grain size below 
100 nm, has been proved to be one of the fundamental options for mitigating these deleterious effects over 
the past 50 years []. The novel radiation resistance of NCs is apparently manifested by the formation of a 
defect denuded zone near the GB, as observed in irradiation experiments of metals [] and oxides []. Quite 
recently, the NCs iron (Fe) with high thermal stability have been prepared successfully [] and see a 
promising application in extreme irradiation environments []. 

Until now, great effort has been made to reveal the atomic origin of the radiation-tolerance of NCs, 
focusing on the fundamental interactions of the V/SIA with the GB in several metals []. It is found that the 
GB has two basic roles either as a trapping center for the SIA/V or as the annihilation catalyst for the SIA-V 
recombination. On the one hand, the trapping of the SIA/V by the GB suppresses the formation of the large 
immobile defect clusters in the grain interior, leading to the low defect density in the grain interior in NCs. 
The mechanism is supported by the rate theory models of the grain size effect on the GB sink strength [], the 
in-situ observation of defect cluster absorption by GBs in Kr ion irradiated NC nickel []. Atomic simulations 
suggest that the low V/SIA formation energy at the GB and low diffusion energy barrier near the GB acts as 
the respective energetic and kinetic driving force for the V/SIA segregation to the GB []. On the other hand, 
several distinct annihilation mechanisms/processes have been proposed in recent years []. In NC copper (Cu), 
the SIA was observed to firstly segregate to the GB and then remit from the GB. Consequently, the V near 


the GB was efficiently annihilated via a low energy barrier process. Such interstitial emission (JE) was 


shown to work near the GB/surface in Fe, tungsten (W), and Cu []. Recently [], it was found that the SIA 
motion along the GB carrying with the low-energy barrier annihilation region around the SIA coupled with 
the diffusion of the V towards the GB could also anneal radiation damage in NC Fe []; the coupling 
mechanism also applies to SIA-V annihilation in W and Cu GBs/surfaces []. The SIA-cluster was found to 
be reflected back into the grain interior as approaching the locally dense GB and annihilate the V near the 
GB in W [], although the SIA-V annihilation could be via the /E and coupling mechanisms near the locally 
loose GB region. 

Despite these advances, significant scientific questions remain unanswered, such as how the 
microstructure near the GB evolves and how the GB relaxes under a harsh cumulative irradiation condition. 
Given that the experimental observation was conducted after a certain radiation dose, the defects 
accumulation and GB relaxation mechanisms are critical to understanding the final experimental 
observations. First, the complex experimental conditions and multiple SIA/V-GB interatomic processes 
render the defects accumulation and GB structural relaxation mechanism hard to access. As illustrated in 
supplementary Fig. S3, the experimental conditions consist of multiple parameters spanning a large range 
(supplementary Table S1), including radiation conditions of temperature (7), radiation dose (rd) and 
radiation dose rate (t), and the grain size (L). The radiation response of the GB is related to the atomic 
processes in the bulk, near the GB and within the GB, and also the process associated with GB motion. In 
different parameter regimes, the system could have various defects accumulation mechanisms, involving 
distinct atomic processes. Second, one also has simulation technique challenge facing the prediction of the 
response of the GB to cumulative irradiation damage of a certain dose due to the possible evolution of the 
SIA/V-GB interaction with radiation dose. The popular simulation methods for exploring cumulative 
irradiation response of NCs are molecular dynamics (MD) and object kinetic Monte Carlo (OKMC) []. In 
MD simulation of cumulative irradiation damage, the repeating ion bombardments with certain kinetic 


energy continuously hit the simulation cell at a certain incident angle at an certain time interval, as 


performed in Fe [], and Cu []. Chen et al. discovered that the SIA-V annihilation at GBs was mediated by 
the formation of chain-like defects []; they concluded that the GBs act as highly efficient defect sinks that 
cannot saturate under extreme radiation conditions. Using similar methods, Jin et al. found the synergistic 
motion of the stacking fault tetrahedra and GBs under repeated bombardments to the GB in Cu []. Our 
simulations on a Fe GB also revealed such mechanism (supplementary Figs. S1 and S2). Although these MD 
simulation results reveal the cascade-GB interaction, the short simulation time interval (an inherent 
limitation of MD) between collision cascades yields a dose rate several orders of magnitude higher than that 
in experiments. As a result, the revealed microstructure relaxation mechanism/process near the GB may not 
occur in the other parameter regime of the real experiments. Alternatively, with the basic SIA/V-GB 
interaction processes as input, the long-time evolution of the radiation-induced microstructure could be 
accessed by OKMC. However, the OKMC model would fail to work if new process emerges as defects 
accumulate with radiation dose. Therefore, in the OKMC model it requires not only the fundamental atomic 
processes to drive the microstructure evolution but also the feedback on the modification of the processes 
from the microstructure subjected to a specific dose of radiation. Such coupling between coarse-grained and 
atomistic simulation is critical for accurately accessing the long-term radiation response of NCs, which may 
also be encountered in stochastic cluster dynamics (SCD) simulation of defects accumulation inside Fe 
grains and in GB []. 

In this work, we proposed a computational predictor-corrector framework for simulating displacement 
damage accumulation and structural relaxation in NC Fe by combining molecular statics (MS), MD and 
OKMC. The method takes advantage of the coarse-grained technique to access the long term evolution of 
the microstructure and the atomic technique to explore the structural transition. The atomic model with 
reproduced defect structure at a radiation dose drawn from OKMC simulations was relaxed by using MD at 
a certain temperature. As a feedback of atomistic simulations to OKMC, the new relaxation process was 


reconsidered in the OKMC model to correct the initial OKMC simulation results. In addition, via calculating 


cumulative distribution function for the defect energy level and its density, the fine and sensitive method 
was developed for characterizing the modification of the whole and local GB structures and the GB's role 
for the V by the SIA segregation. NC Fe was chosen as a prototypical model system since ferritic steels with 
body-centered cubic (bcc) a-Fe as the matrix structure are widely used in nuclear reactors as structure 
materials []. Besides, the basic defect properties, e.g. defect formation energy, diffusion energy barrier and 
defect-GB binding energy in Fe are of the same order of magnitude as that in Cu, Ag and Au []. The present 
results may exhibit universality. In the parameter regime of a high dose rate, a low temperature and a small 
grain size considering experimental conditions [], the SIA overloaded effect was revealed in Fe GBs and an 


alternating evolution of GB structure and the role of trapping and annihilation of Vs was reported. 


2. Computational method 


2.1. Interatomic potential and GB models 
2.1.1. Interatomic potential 

The inter-atomic potential proposed by Mendelev et al. as part of their embedded-atom-method [39] 
(Potential 2 in Ref. [39]) was used to model the inter-atomic interaction in Fe. Tests show that this potential 
correctly predicts the properties of the point defects, e.g., the V/SIA formation energy (1.71/3.53 eV) and 
diffusion energy barrier (0.63/0.33 eV), and the stable configuration of the SIA in the bulk Fe ([110] 


dumbbell) []. The potential has been widely used for studying radiation damage and defect behavior in Fe []. 


2.1.2. Atomic GB models 

The atomic GBs studied in this work include 13 high-angle and small-angle symmetric tilt GBs with 
two kinds of the tilt axes of [0 0 1] and [1 1 0]. These GBs arex 5(210)/[001], x 5(310)/[001], x 
13(510)/[001], = 13(320)[001], x 25(430)[001], x 85(760)[001], = 85(13 1 0)/[001], x 3(111)/[110], 
> 3(112)/[110], x9(114)[110], x9(221)[110], x 11(113)/ [110], x 11(332)/ [110]. The GBs have been 


frequently modeled []. The creation and relaxation procedures are identical to the previously reported 


procedures [29]. Periodic boundary conditions were applied in the two directions parallel to the GB plane. 
The vacuum of 15 A was added in the direction normal to the GB plane and the size in the direction is about 
30 À. The GB system energy was minimized through the rigid-body translations of one grain relative to the 
other, followed by the full atomic relaxation in all three Cartesian directions at 0 K. After relaxation, the GB 
energy density (£,,) per unit area, as defined in Ref. [], was calculated. The relaxed GB structures are 
shown in supplementary Fig. S4. It can be seen, different GBs have distinct structures. Similar GB structures 
are obtained in other work []. It is found that the twin boundary of 23(1 1 2)/[1 1 0] (supplementary Fig. 
S4(1) has the lowest E,, of 0.26 J/m? due to smaller distortion compared to other GBs; similar 
structure-energy relation was found in W []. The GBs of 23(1 1 1)/[1 1 0] (supplementary Fig. S4(h)) and 
29(2 2 1)/[1 1 0] (supplementary Fig. S4(k)) have higher £,, of about 1.3 J/m?. In the main text of the article, 
these GBs will be shown to have rather different defects accumulation and structural relaxation properties, 
which is correlated with E,,. The stability of the GB structure was also investigated by performing MD 
simulations at 300 K for about 50 pico-seconds (ps). The velocity-Verlet algorithm was employed to solve 


the Newton motion equation. The time-step was set 2 femto-seconds (fs). 


2.1.3. Coarse-grained GB models 

For OKMC simulating defects accumulation in NCs, the coarse-grained GB models of the 
square-shaped cubic cells are built (Fig. 1(a)). Periodic boundary conditions were applied in the three 
Cartesian directions. Based on the previous atomic calculations of the basic SIA/V properties near the GB[], 
the whole model was divided into the bulk region, the region nearby the GB and the GB plane. In the three 
regions, the defects have corresponding energetic and kinetic properties, simplified by using different 
physical pictures in the following text. In the OKMC model, the different regimes of the t-L-T space were 
accessed with specific defect-GB interaction parameters. The values of these parameters of t, L, T were 
assigned based on experimental conditions [] in order to make comparison of the theoretical simulation 


results with the experimental results []. Parameter t was assigned to be three levels of 10°, 10? and 10? 


dpa:s !; L was given of 10, 50, 100 and 200 nm; T was given of 100, 200, 300 and 600 K [E38]. The 
purpose of choosing different values of the relevant parameters is to explore parameter-regime dependence 
of the defects accumulation and SIA/V-GB interaction pictures. The other defect-GB interaction parameters 
were obtained from MS calculations. For a large L, it is hard to reach several dpa in a three-dimensional (3D) 
GB model due to the large number of objects included in the model. The two-dimensional (2D) GB models 


were thus also built, while in 3D models L was only given10 nm. 


2.2. Framework of the predictor-corrector method for simulating displacement damage accumulation and 
structural relaxation in NC-Fe 

The proposed predictor-corrector method combines the MS, OKMC and MD, as illustrated in Fig. 1(b). 
The MS was adopted to parameterize a set of unit processes/events. In addition to the defect production, the 
fundamental atomic processes include diffusion of the V,/SIA,, annihilation of the V4;-SIA;,, clustering of 
the V,/SIA,, dissolution of the V,/SIA, in the bulk, near the GB and within the GB, and emission of the 
V,/SIA, from the GB. These basic processes drive the microstructure evolution from one basin to another. 
Since all the unit processes about the behavior of a single V,/SIA, or V4-SIA, interaction near the GB 
known have been incorporated into the model, it is expected that the model is accurate enough to describe 
the structural evolution as the radiation dose is lower than a certain level. The OKMC model at this stage 
gives the SIA,/V, distribution and concentration with radiation dose, as illustrated in Fig. 1(b) by the curve 
of Cs, y~ dose. The OKMC modeling method based on MS-characterized basic processes during this stage 
was termed as the predictor. 

The next step is the corrector. The curve Csj4 y— dose was sampled to obtain the corresponding defects 
distribution at a specific radiation dose and at an elevated temperature. The atomic model with reproduced 
defect structure at a radiation dose was relaxed by using MD at a certain temperature. And then, the defects 
properties of the model were calculated by the MS. After the procedure, the Cs;4 y dependent Egg and 


SIA/V-GB interactions were obtained, as illustrated by the curve of Egg — Csr, y (Fig. 1(b)). Once new atomic 


process was found, the process was again incorporated into the OKMC model to correct the initial OKMC 
prediction of the microstructure evolution based on the new set of unit processes of the V,/SIA,. In principle, 
the prediction-correction procedure is repeated. The details of the predictor-corrector method are described 


as follows. 


2.2.1. MS calculation of the V/SIA properties near the GB 

MS was used to calculate fundamental interaction parameters of the SIA,/V, with the GB, including 
V/SIA formation energy ( E/,,,) [], V/SIA-GB binding energy, V/SIA diffusion energy barriers near the GB 
and within the GB, V-V/SIA-SIA binding energy within the GB. Results for 5(310)/[001] are summarized 
in supplementary Fig. S5. Based on the defects formation energy, the stable configurations of the defects at 
the GB were meanwhile found; the V/SIA configurations at the other several GBs are shown in 
supplementary Fig. S6. Note that, the stability of the site for the V/SIA was also examined at 300 K by using 
MD. It is found that the site only for di-V/SIA obtained at 300 K is slightly different to that at 0 K. The 
energy difference is less than 0.2 eV. The V/SIA segregation energy, characterizing the energetic driving 
force for a V/SIA to segregate to the GB from the bulk, was defined as the reduction in the system energy as 
the defect migrates to the GB from the grain interior, which equals the bulk defect formation energy minus 
the defect formation energy at the GB. In some other references, the segregation energy is also termed as the 
V/SIA-GB binding energy. 

The energy barrier for the V/SIA to overcome to be reemitted from the GB ( £^) was also defined. 
The barrier equals the sum of the V/SIA diffusion energy barrier in the bulk and the V/SIA segregation 
energy. The time needed for the V/SIA emission from the GB (t,,,,,) was calculated via tpp —v,'e^ PT, 
where prefactor v,was given 10'7/s, Boltzmann constant kg has a value of 8.517x10? eV/K, and T is the 
system temperature. To show whether the V/SIA tends to be clustered within the GB, the binding energy of 


the V-V/SIA-SIA at the GB ( £j" 9^) was defined. E/ ^" ^" was calculated by E; ^. =2E} su -Ef isu, . By 


this definition, a positive value of the binding energy means attraction between two Vs/SIAs, while a 


negative value indicates the repulsion between the two defects of the same type. 

The standard nudged-elastic-band (NEB) method [38] was employed to calculate the V/SIA migration 
energy barriers near the GB and along the GB. As calculating the V/SIA migration along the GB, the GB 
was termed as a two-dimensional plane. Two non-equivalent migration paths were determined. To do this, 
the energy barriers along all the possible migration paths connected from the initial site to the other two 
nearest sites were calculated. For example, for the V diffusion from site A to B, the common nearest site C 
of sites A and B was firstly searched. The energy barriers along the paths A-B and A-C-B were then 
calculated; the path with the minimal energy barrier was chosen as the most kinetically favorable one in the 
direction from A to B. Similarly, the migration with low energy barrier on the other direction was calculated. 
The states that have been relaxed for calculating defects formation energy were employed as the possible 
initial and final states of a diffusion process. The dense images of 28 intermediate images were inserted 


between the initial and final states by using the linear interpolation method. 


2.2.2. Model assumption and parameters for the OKMC simulating V/SIA accumulation in NC Fe 
Based on the spherical approximation of an object (Fig. 1(d—e)), the radius of an object of the V,/SIA, 

was defined as follows: 

Am 3 a. 

— R? =n—,r(n)=R+rshell , (1) 

3 2 
where apis the lattice constant of 2.855 A. As shown in Eq. 1, the radius of an object (r(n)) includes two 
parts: the core part R and the shell part rshell. Considering the large stress filed of the SIA, in the bulk, 
rshell for a bulk SIA, was given a larger value than that for the bulk V4. By definition, the radius for a 
SIA,/V, in the bulk and within the GB is given by Eqs. 2—5, as follows: 

3% 3 a 


r(bulk,V ,n) = 5.0 + a, x eo) xn, (2) 


r(bulk, SIA,n) = 7.5+ a, x =" xn, (3) 


r(GB,V,n) =5.0+a, x oy" xn’, (4) 


r(GB, SIA, n) = 5.0 a, «Cm xn? , (5) 

In Eqs. 2-5, the radius is in unit of A. Note that, only when the two defects are both in the bulk or both at the 

GB, the defects were viewed as spherical objects. As we treated the interaction of the V, nearby the GB with 

the SIA, at the GB, the SIA, within the GB was given a limited radius such as 10.0 A because of the 
localization of the SIA composing the SIA, at the GB (Fig. 1(d)). 

In the OKMC model for simulating cumulative irradiation damage in nano-crystal Fe, for a transition 

with an energy barrier of Em, the rate r was calculated viar =v, exp[- E, /(k,7)]. The detailed procedure for 


calculating r of different types of events is as follows. 


(a). Defect production 


For a given dose rate t and grain size L, the SIA/V production rate (rcVI) was calculated by 


atom 


3 
rcVI = Due y , where the volume of the calculation cell (v,,,,) is L^ and the volume of an atom (V, ) is^ 


atom 


2 
Note that, in a 2D OKMC GB model, v,.,,is L’ and y, is^. As the defect production event was chosen 


during the implementation of the OKMC program, a pair of V and SIA was placed randomly in the 


calculation cell. 


(b). Diffusion of the V,/SIA, in the bulk, near the GB and within the GB 

Within the three regions, the V,/SIA, migrates over half the lattice constant each time in one direction. 
Diffusion of the V,/SIA, in the bulk was set to be 3D random walk. Only the small V, with n less than 5 was 
set to be possibly mobile at an elevated temperature (at a high temperature, such simplification may be 
insufficient to correctly predict the microstructure evolution). The energy barrier for the small V,/SIA, 
diffusion in the bulk was from Ref. []. Considering higher mobility of the SIA, than the V,, only SIA, with n 
larger than 10 was set to be immobile, while the barrier for the SIA, (3«n«10) diffusion was set to be that 


for SIA; (the size limit set here for the mobile V,/SIA, ensures the preferential absorption of the SIA, by the 


GB over the Vn). The diffusion energy barriers for the bulk SIA,/V, are summarized in supplementary Table 
S2. 

The physical picture for the V,/SIA, behavior near the GB considered in the OKMC model includes 
GB-enhanced diffusion of the V,/SIA, nearby. Diffusion of a V,/SIA, near the GB was considered to be 
one-dimensional walk. Only the V, with n less than 5 was considered to be possibly trapped by the GB (such 
approximation may not hold true at a high temperature, as proofed by the absorption of defects clusters by 
the GB in experiments []). The energy barrier for the GB trapping of the V, was set to be a half of the barrier 
for the V, diffusion in the bulk given the accelerated diffusion of the V, by the GB []. Since the SIA, near 
the GB was spontaneously trapped by the GB, the energy barrier for the SIA, migration near the GB was set 
to be zero. The energy barriers for the SIA,/V, migration near the GB are shown in supplementary Table S2. 
The up limit of n was set 500. The trapping radius of the GB for the V, and SIA, was set to be 5 and 10 A, 
respectively. To execute the trapping of a V,/SIA, by the GB, the V,/SIA, was put at the GB plane. In the 
present square-shaped model, there are three GB planes. The plane with the shortest distance from the 
VSIA, traps the V,/SIA,. The simplified picture is illustrated in Fig. 1(d) and (f). Diffusion of the V,/SIA, 
within the GB was modeled to be 2D; only the single V/SIA was considered to be mobile within the GB (Fig. 
1(e)). The energy barriers for the V,/SIA, migration on the two nonequivalent directions within the GB are 


shown in supplementary Tables S3 and S4. 


(c). Annihilation of the V4;-SIA;? in the bulk, near the GB and within the GB 

As both the Vn; and SIA, are located in the bulk or within the GB, the energy barrier for the 
annihilation of the V5; with SIA;? was set to be zero. For the annihilation of the V4; near the GB with the 
SIA,2 at the GB (Fig. 1(f)), the considered annihilation mechanisms include the annihilation induced by the 
SIA-emission [] and the V-hop []. Given the results of MS calculations that reveal the formation of the 
spontaneous annihilation region surrounding the SIA, at the GB [], the radius of such region was set to be 5 


A. Out of the region, the energy barrier for the recombination varies lineally with the distance of the V; 


from the border of the annihilation region; the barrier converges to the value for the V4; migration in the 
bulk. The approximation picture is illustrated in Fig. 1(f). After annihilation, the defect with a larger size 


survives with the reduced size of|n- 2|. 


(d). Clustering of the V,/SIA, in the bulk and within the GB 

The energy barrier for the V clustering in the bulk was set to be half that for the single V migration in 
the bulk. Only the cluster V, with n less than 5 was allowed to be trapped by another V4. Considering the 
high mobility of the SIA, the energy barrier for the SIA clustering in the bulk was set to be zero. 

For the SIA/V clustering within the GB, the energy barrier was set to be the half of the minimal barrier 
for the V/SIA migration along the GB if the binding energy of the V-V/SIA-SIA is positive; otherwise, the 
barrier for the clustering was assigned to be the V/SIA migration energy barrier along the GB minus the 
V-V/SIA-SIA binding energy. Only the single V/SIA was allowed to be trapped by another cluster. The 


binding energies are given in supplementary Tables S5 and S6. 


(e). Dissolution of the V,/SIA, in the bulk and within the GB 

Each time, only one V/SIA was dissolved from the V,/SIA,. The energy barrier for the V/SIA 
dissolution from a V,/SIA, in the bulk and within the GB was set to be the sum of the migration energy 
barrier for the single V/SIA and the binding energy of the V-V,/SIA-SIA, in the corresponding region. Only 
small clusters with n less than 4 were allowed to be dissolved. To execute the V/SIA dissolution from a 
V,/SIAy, a V/SIA composing the V,/SIA, was put at a site that is d away from the V,/SIA,; distance d is the 
sum of radius of the Vn-;/SIA,1 and that for the V/SIA plus 10 A to avoid the re-trapping of the V/SIA by 


the V,/SIA, in the next OKMC step. 


(f). Emission of the V,/SIA, from the GB 
The energy barrier for the V,/SIA, emission from the GB is the sum of the V,/SIA, segregation energy 


and the V,/SIA, migration energy barrier in the bulk, as listed in supplementary Tables S7 and S8. Each 


time, only one of the V/SIA of the V,/SIA, was emitted from the GB. During executing the emission event, 
the emitted V/SIA was put the bulk site that is d away from the GB. To avoid the subsequent frequent 
re-trapping of the V/SIA by the GB, the distance d was set to be the GB trapping radius for the V/SIA plus 5 
À. Note that the emission described here is different to the SIA-emission that induces the annihilation of the 
V nearby [], which is the atomic detail for the V-SIA annihilation near the GB. 

We found that, as all the atomic processes in the bulk, near the GB and along the GB were incorporated 
into the model, the long timescale or high radiation dose could not be accessed in some cases. The 
preferentially segregated SIA at the GB was observed to frequently migrate within the GB at high 
temperature, and consequently the V in the bulk region nearly remained immobile. In these cases, the 
diffusion of the V/SIA along the GB was inhibited. In the one hand, such picture corresponds to the GB with 
large V/SIA diffusion energy barrier along the GB. In the other hand, in the case of the very small barrier for 
the V/SIA along the GB, such assumption implies that the SIA quickly adjusted its position during the 


occurrence of other events. 


2.2.3. Details for MD simulating GB thermodynamics under different levels of V/SIA concentration 

With the defects distribution in hand, the structure and energetics of the GB under different levels of 
V/SIA concentration were investigated. After loading a certain number of the V/SIA into the GB, the system 
was firstly relaxed at 0 K by using the steepest descent method, followed by the MD relaxation at 300 K. 
After the relaxation, the E..,, and corresponding structure with a certain V/SIA concentration were obtained. 


To obtain the statistical results of the GB structure and E, , 


10 distinct MD simulations were performed. 

To access more discrete V/SIA concentration values at the GB, the number of the structural unit in the 
two directions parallel to the GB varies. Depending on the size of the minimal structural unit of a GB, the 
number lies in the range 2—5. In addition, different SIA loaded styles were considered given the short time 


scale of MD and the resulted possibly metastable state of the GB with the loaded defects. In Stylel (Fig. 


1(g)), the site for the single SIA/V occupation within the GB is first found based on the calculation of the 


energetic property of the SIA/V near the GB. Then, a certain number of SIAs/Vs are put at these sites to 
produce the required SIA/V concentration. It could expect that such style would generate a defect structure 
close to the stable state under a low dose. In Style2 (Fig. 1(h)), the SIA/V is put near the GB one at a time, 
followed by the MD relaxation at a certain temperature, e.g. 300 K. The sequentially added SIAs/Vs are 
finally trapped at the GB. 

To characterize the GB structure with different concentrations of defects, some other more effective 
methods were developed in addition to the Egg and the visualization method. Because, the GB was found to 
develop to many rather different local structures under accumulated radiation damage, while Egg measuring 
the stability of the whole GB structure, is insensitive to the local GB structure. Based on the V formation 
energy near the GB with a certain SIA concentration at the GB (C77), we calculated the cumulative 
distribution function of the defects energy level and the density of a specific energy level. The distribution 
function and density were found to be sensitive to the distinct local GB structures. 

For a given GB with a certain concentration of SIAs, the V formation energy at the sites within the 
range with a radius of 10 À was calculated. Results suggest that, such width of the range is large enough for 
covering the region affected by the GB. The defects energy level cumulative distribution function ( C(E/) ) is 


calculated as 


YN < ED) 
CE —— , (6) 


tot 


where N; is the number of sites with E/ less than E/ , and N; is the total number of sites in the calculation 
region. C( E.) could give the cumulative probability of the V with a certain energy level ( E/ ). Generally the 


V near the GB has multiple discrete levels of the formation energy near the GB [], which makes the 


occurrence of the terrace on the C( £7) curve. The £/ of the end point of the terrace defines the V energy 


level. Note that, since only the sites near the GB are underlined, the sites in Eq. 6 are restricted to ones with 


Ej, less than the bulk V formation energy a certain value, e.g. 0.3 eV. Otherwise, the bulk sites take a too 


large fraction of the calculated region to effectively analyze the evolution of the site property (trapping or 


annihilation) with radiation dose. 

According to the physical meaning of the E/ (e.g. a negative value of E/ often indicates the instability of 
the site and the spontaneous annihilation of the V with the SIA []), the C( £7) was further divided into 
several regimes. The site density within a specific regime (C,,, ) was calculated as 

Cs = Nux[ICG?) - CCZ))/ Sas , (7) 


where E; and E/ are the lower and upper limits of the E/ in the regime, and S,,is the GB area. The density 


could give the information on the relative importance of a certain atomic process. 


3. Results 


3.1. OKMC simulations of microstructure evolution near 35(3 1 0)/[0 0 1]under cumulative irradiation 
With the fundamental interaction parameters of the SIA/V with the GB of X5(3 1 0)/[0 0 1] (as 
summarized in supplementary Fig. S5) in hand, OKMC simulations were firstly performed to explore 
defects behavior in the GB system under cumulative irradiation. The parameters L, T and t are varied. To 
characterize the evolution in the microstructure driven by the defects accumulation in the system, we 
calculated the defects concentration and distribution near the GB and at the GB, as well as their variation 
with radiation dose. For the L of 10 nm, results are mainly presented in Fig. 2, while results for other cases 
are shown in supplementary Figs. S7 and S8. We found that defects accumulation behavior depends on the 
parameter regimes of L-7-: space, as illustrated in Fig. 2(j) for the defects concentration level in different 
parameters regimes. Apparently, in most of cases the SIA concentration in the bulk ( C7^ ) is lower than the 
bulk V concentration (C/"); it is in agreement with the atomic calculation predicting the preferential 
segregation of the SIA to the GB over the V []. Particularly, in a system with a small L, e.g. 10 nm, the 
extremely small value of Ci" was observed, indicating that no SIAs survive in the bulk under such condition. 
It may be due to the too small value of L for the SIA to aggregate to large stable SIA,s. Such phenomenon 


weakly depends on T and t possibly because of the very low diffusion energy barrier of the SIA near the GB. 


In contrast, T and c significantly affects V accumulation near the GB. At a high T and/or low r, the system 
has a low C7" . Especially, as L is smaller than 100 nm, T is higher than a certain level e.g. 300 K and 1 is 
lower than a level e.g. 10? dpa-s`', there is nearly no V surviving near the GB, as suggested by the quite low 
level of C" , This is in agreement with the expectation that the GB system with a small L should well 
self-heal the radiation damage subjected to under a certain kinetic condition. It could be attributed to the 
enhanced diffusion/annihilation of the V [] and reduced formation probability of immobile V, near the GB 
in the corresponding parameter regime. It also implies that the radiation resistance of NC Fe in the regime 
could be explained based on the fundamental SIA/V segregation and annihilation processes included in the 
OKMC model (Fig. 1(b)). 

Besides these observations consistent with our intuition, we however also observed exceptional results. 
In a certain regime of the L-T-t e.g. a low T and high t, C” increases with radiation dose, even if when L is 
as small as 10 nm (Fig. 2(a—c)); it implies that the V accumulates in the bulk with radiation in the considered 
parameter regime. Meanwhile, SIAs accumulate at the GB, as indicated by the increase in C3 with dose (Fig. 
2(e—g)). Besides, further examination of the defect distribution near the GB (Fig. 2(k-n) and supplementary 
Fig. S8) suggests that after a certain dose, the C/" reaches a peak in the bulk region quite near the region 
with a very low C7" in the vicinity of the GB; C7" also decreases with the increase in the V distance to the 
GB due to the periodical property of the OKMC model. These results suggest that the OKMC model 
containing annihilation and segregation processes of the SIA/V predicts the accumulation of radiation 
damage in the bulk region near the GB and at the GB, which seems to contract the self-healing property of 
the small-sized GB system with low-energy barriers for the V/SIA diffusion and annihilation nearby []. Then, 
we are naturally faced with the following issues: whether both C?" and C would continue to increase as 
radiation proceeds and how the accumulated defects modify the GB structure and consequently affect the 
SIA/V behavior near the damaged GB (Fig. 2(k-n)). If the trend of C7" and C27 fails to increase after a 


critical dose, it requires the clarification of new atomic processes or mechanisms occurring, which have to 


be incorporated into the OKMC model to correctly predict the radiation damage accumulation in NC Fe 


together with the ones in Fig. 1(b). 


3.2. MD simulations of the structural relaxation of the GB 35(3 1 0)/[0 0 1] with accumulated defects 

To answer the questions raised above, we sampled the OKMC-calculated C7,,/ C -dose curve (Fig. 
4(a-d)) and revisited to the defects accumulation and structural relaxation mechanisms near the GB on the 
atomic scale by using MD. We concentrated on the regime of the L-T-t parameter space where the radiation 
damage could not be well healed as predicted by initial OKMC simulation results (Fig. 2). Although the 
OKMC simulations reveal that Vs accumulate near the GB and SIAs aggregate at the GB, firstly only SIAs 
were contained in the atomic models at the GB to explore how the GB relaxes after it absorbs defects. To do 
this, a certain number of SIAs were loaded into the GB by using the method as described in section 2.2.3; 
the relation of the C2, to radiation dose and L-T-r could be obtained by referring to A/ in Fig. 2(e). The GB 
with SIAs was relaxed sequentially at 0 and 300 K. The Egg was then calculated. By analyzing the change in 
Egg and the corresponding GB structural variation process, we revealed the GB structural response to 
cumulative irradiation. 

Figure 3(a) shows the Egg profile as a function of cS. The profile could be divided into three stages 
with different features of the E,, changes with CZ/. From Co to about C., the Egg nearly increases linearly 
with C7, , as reported in Ref. [] for one Fe GB of X17(5 3 0)/[0 0 1]. Such tendency could be derived from 
the definitions of the Egg, E7" and E; (supplementary text S1). For the isolated SIAs at the GB in the case 
of a low C; produced by a low dose irradiation, the slope of the curve only depends on the £?“ at the GB. 
Yet, with more SIAs introduced into the GB, C$ is over than the level of C$ (Fig. 3(a)), the growth in Egg 
slows down, as shown in Fig. 3(a) from C; to Ce. During this growth stage from Cy to Ce, a noticeable 
change in E,, was observed. The increase in Egg due to the SIA accumulation at the GB could be readily 
understood based on the contribution of the defect formation to Egg. However, as CH increases to a certain 


level ( C27; in Fig. 3(a)), the Egg no longer increases but exhibits a drop trend with the increasing in the C27. 


The decrease trend holds from Ce to Ce. The different dependence of Egg on the C, from C; to C, to that 
from C, to C; was termed as the SIA overloaded effect. Of course, the overloaded SIAs at the GB would 
also affect the GB structural transition and the subsequent defect-GB interactions. After C,, the Egg seems to 
increase and decrease periodically (*Cyclel" and “Cycle2” in Fig. 3(a)), although it changes smaller than 
that from C; to C, (*Cycle0" in Fig. 3(a)). 

To reveal the GB structural evolution and corresponding relaxation mechanism during the different 
E, ~ Cg, stages in Fig. 3(a), we examined the relaxed GB structure after each loading of the SIAs into the 
GB with focus on the effect of the overloaded SIAs on the GB structure. Several typical GB structures with 
different levels of C$’ are shown in Fig. 3(b-j). We found that, from C, to Ce, the GB containing a low 
concentration of SIAs remains its basic structure (Fig. 3(b) and (c)) as a pure one without defects 
(supplementary Fig. S4(b)); in other words, the structural unit along the GB is not modified by the absorbed 
SIAs. Meanwhile, the loaded SIAs also maintain their identity at the GB, as shown in Fig. 3(b) and (c) by 
the red spheres with high potential energies. The structural observation here is consistent with the linear 
dependence of Egg on C% (from Co to Ce in Fig. 3(a)), as predicted by the supplementary text S1 on the 
theoretical derivation of the relation of the E,, to C7, . It also implies that, in this SIA concentration range or 
corresponding dose range, the OKMC modeling results (Fig. 2) are correct due to the fidelity of the 
important SIA-induced annihilation process incorporated into the OKMC model. 

Nevertheless, from C; to Ce, parts of the loaded SIAs were found to disappear after the GB relaxation. 
Particularly, the new GB structure was formed as more SIAs were introduced into the GB (Fig. 3(e)). The 
structure with the Egg of 1.35 J-m? is obviously different to that for the pure GB with the Egg of 0.98 J m°’. 
The atom separation distance in the new GB system also exhibits unique quasi-continuous distribution 
compared to the discrete distribution in a pure GB (inset figure in Fig. 3(e)). In the latter section, we will 
show the new GB structure has its special defects properties. In some cases of MD simulations, we also 


observed the formation of the other new GB structures induced by the introduced SIAs (supplementary Fig. 


S9(a) and (b)). Similar formation of the new GB phases induced by the defects accumulation at the GB has 
been reported in several face-centered cubic (fcc) GB systems []. Therefore, the interaction among excess 
SIAs and the new phase formation could give rise to the reduced growth rate in Egg from C; to C, (Fig. 3(a)); 
the attraction between SIAs and the transition of the loaded SIAs to components of the GB release the 
system energy. From C, to C,, the GB gradually recovers to its original structure, accompanied by the 
migration of the local GB parts (Fig. 3(f)). More specifically, some Vs in the form of V-SIA-V complexes 
appear near the sites of the lost SIAs at the GB. Finally, the SIAs and Vs coexist at the GB. The recovery of 
the GB structure via the GB migration and V formation could thus lead to the decrease trend of the Egg with 
C$, from C, to C, in Fig. 3(a). The Cf -dependent GB structural relaxation behavior suggests that it requires 
different OKMC modeling schemes before and after the critical radiation dose. For example, as CH 72.5 
nm? (Fig. 3(a)), we could obtain the critical dose of about 0.75 dpa for a given radiation condition of L=10 
nm, 7-300 K, and 7-10? dpa:s'! (Fig. 2(g)). 

Since from Co to Cg, the SIA accumulates at the GB, followed by the GB relaxation and recovery, the 
period from Co to C; could considered as a cycle, as marked by Cyc/e0 in Fig. 3(a). As the GB in Fig. 3(g) 
containing a certain concentration of Vs and SIAs is selected as the starting structure for absorbing the 
subsequent SIAs, it could expect the similar behavior of the GB to that in Cyc/e0. For instance, we observed 
the formation of new GB phases, although it is not so frequent due to the existence of the Vs. We more 
frequently observed SIA-V transition during the GB relaxation and also found the enhanced GB mobility. 
Therefore, as the GB continuously absorbs SIAs, new atomic process occurs; the overloaded SIAs could not 


only affect Egg variation trend but also modify the GB structure. 


3.3. SIA-V transition process near 35(3 1 0)/[0 0 1] 
To reveal the transition process for the SIA to the V and the underlying physical reason, we analyzed 
the atom displacements field near the GB during the structural relaxation. Fig. 4 shows the atom 


displacements of the GB of Y5(3 1 0)/[0 0 1] with a certain concentration of the SIA at 300 K for 50 ps 


relative to the initial GB state. As the SIA was put near the GB in Style2 (Fig. 1(h)), the atom displacement 
was induced by both the SIA segregation and GB relaxation. Furthermore, the atom displacement due to the 
formation of new GB phases (Fig. 3(e)) relative to a pure GB is also involved in the transition. The atom that 
are aligned in one column near the GB in the final state are however not well aligned in the initial state, e.g. 
atoms in columns C and D (Fig. 4(b)). Consequently, the displacement field is complex (Fig. 4(a)). Yet, we 
still could observe that it lacks one atom in column G/ in the final state, forming a V at the GB; meanwhile, 
in the initial state the atoms or SIAs are not fully occupied in one column (Fig. 4(c)). We speculated that the 
overloaded SIAs at the GB occupied one atom column after transition, forming regular atom columns but 
with a certain number of Vs due to the incomplete occupation of the atoms in the column. 

Such picture was further validated in another simple model where SIAs were loaded into the GB via 
Stylel (Fig. 1(g)); the SIAs were intentionally put at the SIA site at the GB to avoid the coupling of the 
displacement fields of the SIA segregation and that for the SIA-V transition. After comparing the GB 
position before and after relaxation, we found that the GB migrated to the nearest state (Fig. 4(d)). Besides, 
the displacement fields near the GB core are significantly different from that far from the GB. Some regular 
displacement patterns correspond to the intrinsic motion of the atoms during the GB migration, while others 
were induced by the transition from the SIA to V. Fig. 4(e) shows the corresponding atoms in the initial and 
final states. The above transition picture was confirmed, as illustrated in Fig. 4(f). The SIA; in a GB structure 
unit was initially aligned along the tilt axis. The number of defects locations along the tilt axis in a 
computational cell is N. As the number of SIAs i approaches to N, the GB structure was transformed to the 
structure with a V; (j=N-i) after tens of ps at 300 K. Meanwhile, the j atoms on the opposite GB side move to 
the GB plane and become the SIA; leaving a Vj behind. The structure of a V;-SIA;-V; forms (the MS 
calculation suggests that the energy of the GB with a V-SIA-V complex is about 0.1 eV lower than that for 
the GB with a V). As for the transition from a new GB phase (Fig. 3(e)), it could be considered as transition 


of the atoms in Fig. 4(d) and (e) with extra displacements related to the new GB phase. 


The energetic driving force for the SIA-V transition was then analyzed. In supplementary text S1, we 
have derived the relation of the GB energy variation (A£,, ) to the defects formation energy ( £5, ) and the 
defects concentration (C ) at the GB approximately as AE, = E94 xC. Therefore, the two factors of E94 and 
C influence Az,, and further affect the transition process. As shown in supplementary Fig. S5, Ey, is lower 
than E; , implying that it is energetically unfavorable for the SIA to transit to the V. Only when the excess 


SIAs are transited to the minority of Vs could Az,, be greatly reduced during the transition. 


3.4. Evolution in the GB’s role for Vs under cumulative damage 

OKMC simulations in Fig. 2 and supplementary Figs. S7 and S8 suggest that at a specific radiation 
condition the V mainly survives near the GB. Meanwhile, Figs. 3 and 4 show the variation in Egg and GB 
structural evolution during the cumulative irradiation. Thus, it is of meaning to explore the corresponding 
variation of the GB's role as a sink for the V during the cumulative radiation. Previously [], it has revealed 
that the GB mainly has two kinds of roles for the V either as the trapping center for the V or as the catalyst 
for the V-SIA annihilation, closely related to the sign and magnitude of E/. To reveal the evolution in the 
GB's role for Vs under accumulated radiation damage, we calculated the E/ near the GB with different 
concentration of SIAs as well as the cumulative distribution function of the E/ (C(E/) ). As shown in Fig. 
5(a) the GBs with several levels of C% (atomic structures in Fig. 3(b-j)) have characteristic C(E/) . 
Depending on the level of the C2, , the C(E/) for these GBs could be divided into three categories. The pure 
GB (supplementary Fig. S4(b)) only has £; distribution in the regime of the large positive E; ([ Ej , E^ |). 
The GB in Fig. 3(e) has a quasi-continuous distribution of the E/ , starting from zero point to the £7" ((0, 
Ej" |). The GBs with the other levels of theC;; have similarC(E/), featured by the high contribution of 
large negative E/. Actually, the different features of theC(E/) suggest the changes brought about by the 
accumulated SIAs at the GB related to the V-GB interaction in terms of the V-GB interaction nature and 
binding strength. Specifically, the sign of the £/ is positive for the V near a pure GB, and then the sigh 


becomes negative as for the GB loaded with a certain number of SIAs. However, above a critical value of 


the C£ , the sign for most of sites is positive again. We also found that the V-GB binding strength is greatly 
enhanced, as manifested by the occurrence of small but positive E/onC(E/). This could arise from the 
enhanced effect in the new GB phase with a high atomic density and the binding of the V to the V transited 
from the SIA. The variation in the V-GB binding energy will affect the basic behavior of the V near the GB. 
Therefore, as the C$ varies, the V-GB interaction nature would be modified; meanwhile, the GB's role for 
the V may dynamically change. 

To better characterize evolution of the V-GB interaction, we further categorize the sites near the GB 
with different C of SIAs based on the feature of the C(E/). Considering the physical meaning of E/, the 
C(E/ ) was further divided into three groups (Reg/, Reg2, and Reg3 in Fig. 5). The Vs in each regime has a 
Ej within a certain range, which was considered as an general energy level for the V. In Reg/, the V has 
clear energy levels of large positive value and small displacement after relaxation of the GB with the V (Fig. 
5(b)). It implies that after the bulk V enters Reg/, the system energy reduces but the V remains locally at the 
site. Therefore, the region could be considered as a trapping region for the bulk V. In Reg3, the V also has a 
large magnitude of E/, but the value is negative; the atom has a large displacement after relaxation of the 
system with a V in the region (Fig. 5(b)). We have shown that it forms the spontaneous annihilation region 
near the trapped SIAs with such E feature []. And thus, Reg3 could be defined as the recombination region. 
We also found additional feature of E/ that has a very small value close to zero in Reg2, e.g. the results for 
the C2, of C, in Fig. 5(a). It may be related to the locally disordered structure (Fig. 3(e)). 

To further qualitatively describe the probabilistic V-GB interaction and its evolution with the C2; , we 
further introduced the concept of the density of the V energy level (characterized by E/) in the three regimes. 
Combining the physical meaning of each regimes and the density variation with the C7 , we could obtain the 


evolution information of the GB's role for the V. Fig. 5(c)) presents the calculated density ( C, ) for 


eg1,Reg2,Reg3 
each GB system loaded with different C; of SIAs. As can be seen from Fig. 5(C), Cregi resz regs Varies 


remarkably with C$} , suggesting that C, is more sensitive than E,, (Fig. 3(a)) to the microstructure and 


eg1,Reg2,Reg3 


may serve as a more effect method for characterizing the GB structural evolution during irradiation. 


However, the accumulated SIAs at the GB have totally different effects on Cr and Cres given their individual 


egl 


trends (Fig. 5(c)). Specifically, from J to P/,C 


res increases with Cy, implying a more important role of the 
GB in annihilating the V nearby as more SIAs accumulate at the GB due to the more contribution of the 
atomic processes related to annihilation e.g. JE to alleviate the Vs near the GB with the SIAs aggregated at 
the GB. Such enhanced effect of the GB's role in annihilating defects could be readily understood based on 
the natural extension of the annihilation region near the multiple SIAs compared with the single SIA. 
Meanwhile, the reverse trend for Cp with C7; was observed within the same C; range, e.g. from /0 to S1 
within the first Cycle of Cycle0. It indicates that the GB's role as a trapping center for the V becomes 
weakened. Exactly, parts of the trapping sites for the V near a pure GB were changed to the annihilation sites 
near the GB with a low concentration of SIAs (Fig. 3(a—c)). The reason could be the instability of more sites 
near the GB induced by the more SIAs at the GB. Note that, for a smaller E; (as described in the caption of 


Fig. 5), C,,, would increase with Cy, suggesting that the SIAs at the GB only slightly affect the stability of 


the sites with positive E/. From P/ to B1,C 


regs Feduces greatly with Cy; , indicating the weak role of the GB in 
annihilating the V nearby. This is consistent with the structural observation in Fig. 3(b-e) that at 
BIcorresponding to the Cy; of Ce the loaded SIAs become components of the new GB phase. After S1, Cre 
increases greatly and reaches its peak value at 77, while C,,,has a quite small value. It suggests that near 
the new GB phase, the V interacts with the GB mainly via the trapping process. The extended strain field 
near the GB (Fig. 3(e)) by the large atom density could account for the enhanced trapping for the V. The 
similar trends of C,,, and Cre, were observed from T7 to S2 or from B/ to P2 to that from J to PJ. The 
structural reason for the trends is also similar. In the two cases, the accumulated SIAs at the GB induced the 
transition of the trapping site to the annihilation sites due to the lost identity of the SIAs at the GB. Only the 


initial GBs at 70 (supplementary Fig. S4(b)) and T7 (Fig. 3(e)) have different phases. The special structural 


reason during this concentration range includes the recovery of the GB from T7 to S2. Therefore, within the 


first cycle of Cycle0, we revealed an alternating evolution of the GB's role of trapping and annihilation of Vs 
compared to the previously known unique role of the GB either as a trapping center for the V or as the 
annihilation catalyst []. Furthermore, the alternating pattern or mechanism leads to the different dominated 
roles that the trapping and annihilation play in the V-GB interaction within different C7; ranges or 


equivalently under different radiation conditions. For instance, by comparing the magnitude of Cre and C,.,;, 


we found that the trapping plays a dominated role in removing the V near the GB from 70 to JO, while from 
J0 to P1 the annihilation dominates. The present results suggest that, the previous reported self-healing 
mechanisms for radiation damage in NCs based on the SIA behavior near the GB (e.g. the annihilation 
induced by the JE [] and the coupling of the V migration towards the GB with the SIA moving along the GB) 
may fail under extreme irradiation conditions due to either the new phase formation of the GB or the GB 
structural recovery. 


After Cg, although no fine features were observed onC, compared with that within Cycle0, some 


eg1,Reg2, Reg3 


similar trends were still noticed. The C,.,, for a GB with a certain concentration of SIAs is much higher than 


C,,,, and also than Cp for a pure GB without SIAs, indicating that the trapping is a more basic process than 


annihilation in the GB system exposed to continuous irradiation. It could arise from the formation of the V 
after GB structural transition (e.g. Fig. 3(g)) Basically, C,,, and Cp, show a reverse trend with C5; . 
Exception was encountered during a certain C, range (e.g. from about 17 to 21 nm?) where Cree and C,,,, 


increase simultaneously due to the formation of the V. Finally, some special features were observed on the 


curve of the Cp. with Cz, , which reflects the structural meaning of C,,,.and its connection to GB structures. 


Reg2 eg2 


For example, C, reaches a peak at U1 in Fig. 5(c) where the new GB phase also forms (Fig. 3(e)). Besides, 


eg2 


Cree has a larger value within Cycle] and Cycle2 where the GB consists of V-SIA complexes and is not so 


regular compared to a pure GB. 


3.5. Evolution in the kinetic property of Vs under accumulated damage 


We next examined the V kinetics near several typical GB structures as observed in Fig. 3(b-), including 


a pure one (Fig. 6(a)), the GB with multiple SIAs with well-defined identity (Fig. 6(b)), the GB with a new 
phase (Fig. 6(c)), and the GB after recovery (Fig. 6(d)). We calculated the V migration energy barriers along 
one typical paths near these GBs by using the NEB method (section 2.2.1) []. As seen in Fig. 6(e), all the 
profiles show a decreasing trend with the V closing to the GB, implying the enhanced kinetics of the V near 
the GBs. The difference is also significant. Comparing the results for the GB in Fig. 6(a) and Fig. 6(b), we 
found that the V is finally trapped at the GB along the diffusion and trapping path of path/, while along 
path2 the V would be annihilated at site E2 within the annihilation region Reg3. Furthermore, the energy 
barriers along path2 are smaller than that along path/, indicating that SIAs at the GB promotes the V 
diffusion and annihilation near the GB via reducing the energy barriers. These results are consistent with 
previous results []. Comparing the results along path] and path3, we found that the energy barriers for the V 
diffusion along path3 are also smaller than that along path1, implying that the kinetics for the V near the 
new GB phase is further enhanced than that near a pure GB. Since both sites G/ and F3 locate within Reg/, 
path1 and path3 are diffusion-trapping paths for the V. We also found the promoted diffusion of the V along 
path4 near a GB containing the V after recovery compared to that along path/. Therefore, no matter how the 
GB relaxes after trapping a concentration of SIAs from the bulk region nearby, the diffusion of the V near 
these GBs is always accelerated due to three factors of the accumulation of a specific concentration of SIAs, 
or the new phase formation of the GB or V formation after the GB structural recovery. 

These enhanced processes for the V predicted by the above MS calculations of the V kinetics near the 
GBs were validated by MD simulations of the V behavior near the three typical GB structures at 300 K (Fig. 
6(f-h)). A certain concentration of Vs were initially put near the GB. We found that after 50 ps at 300 K, the 
Vs near the pure GB remain immobile (Fig. 6(f)), while the Vs near the new GB phase were trapped by the 
GB (Fig. 6(g)). The Vs closer to the GB after the GB's recovery also segregated to the GB (Fig. 6(h)). 
Therefore, it becomes easier for the Vs near the GB to enter the GB with a new phase or containing Vs after 


recovery. 


3.6. Defects accumulation and structural relaxation mechanisms in the other GBs 

The above OKMC, MS and MD results (Figs. 1—6) are obtained in the GB of X5(3 1 0)/[0 0 1]. To show 
whether the defects accumulation features and GB structural relaxation mechanisms apply to the other GBs, 
we also performed extended similar simulations in several other GBs (the relaxed pristine structures of these 
GBs are given in supplementary Fig. S4). We compared the atomic and coarse-grained defects behavior in 
these GBs. We found that both the V and SIA distributes within several layers near the GB (supplementary 
Fig. S6). The V locates at the stable site of the GB as a point (Fig. S6). However, the SIA configurations are 
complex (supplementary Fig. S6). In addition to the point configuration at the GBs of X5(2 1 0)/[0 0 1], 
X13(5 1 0)/[0 0 1], £13(3 2 0)/[0 0 1], X85(13 1 0)/[0 O 1], £3(1 1 1)/[1 1 0], X9(2 2 1)/[1 1 0], other 
configurations were also found, including the dumbbell at the GBs of X25(4 3 0)/[0 0 1], 285(7 6 0)/[0 0 1], 
X9(11 4y[1 1 0], X11(1 1 3y[1 1 0], X11(3 3 2y[1 1 0] and the crowdion at the GB of 3(1 1 2)/[1 1 0]. The 
different SIA configurations at the GB may be related to the atom density near the GB core. For example, 
the twin boundary of 23(1 1 2y[1 1 0] with the SIA of a crowdion has the largest atom density near the GB 
core. 

The OKMC parameters also depend on the GB structures (supplementary Tables S3, S4, and S6-S8). 
Based on the MS parameters, particularly the emission energy barrier for the V/SIA from the GB ( E^ / EF) 
and the corresponding time (1, /t,.), prejudge the V/SIA behavior near the GB. We found that the twin 
boundary has the lowest EF . Consequently, only at the twin boundary is ton the order of one hour at 
room temperature, while in the other GBs tis far larger than one hour, particularly for the SIA). Thus, 
the V/SIA could be reemitted from the twin GB once segregated to the GB at an elevated temperature. The 
further OKMC simulations of a single SIA/V behavior suggest that at 300 K, the SIA near all the GBs first 
randomly walks in the bulk region and finally gets trapped at the GB; it is tightly bound at the GB after 
trapping (supplementary Fig. S10(a)). Yet, the V was observed to be occasionally reemitted from the twin 


GB at 300 K after segregation to the GB (supplementary Fig. S10(b)). As the 7' increases to 400 K, both the 


SIA and V were frequently remitted from the GB (supplementary Fig. S10(c) and (d)), particularly for the V 
(supplementary Fig. S10(d)). It could expect that the difference in the defect accumulation rates in the two 
GBs is not significant at a low T. The emission in the twin boundary enhances bulk V accumulation rate and 
decreases the healing efficiency for radiation damage in the bulk region by the GB at the parameter regime 
of a high T and a low t. Therefore, in the parameter regime we are interested in e.g. low T and high z, the 
reemission of the SIA/V from the GB may not affect the SIA/V accumulation and GB relaxation mechanism 
significantly; of course, the critical dose for the new GB relaxation mechanism to occur (point A on the 
curve of Cs; y-— dose in Fig. 1(b)) would be modified once the remission was activated. 

We next mainly explored the curve of Egg~ Csi4y in the other GBs by performing similar relaxation 
simulations to that in Z5(3 1 0)/[0 0 1]. Yet, MD simulations at 300 K suggest that different GBs have 
distinct intrinsic mobility, as measured by the GB migration distance within a time interval (supplementary 
Fig. S11). Some of the GBs remain immobile at room temperature, e.g., the GBs of X5(2 1 0)/[0 0 1] and 
X5(3 1 0)/[0 0 1]. However, the GBs of 225(4 3 0)/[0 0 1], X85(7 6 0)/[0 0 1], 2285(13 1 0y[0 0 1], X3(1 1 
1y[1 1 0], £9(2 2 1y[1 1 0] and X11(3 3 2)/[1 1 0] were observed to move and even disappeared via 
interacting with the free surface in the GB system . We thus only present results for the GB that did not 
disappear after MD simulations of the GB relaxation for 50 ps. Figure 7(a—1) presents the E, profiles as a 
function of C, in these relatively stable GBs. First, it was noticed that the magnitude of AZ,, varies greatly 
in different GBs. Since £,, measures the thermodynamic stability of a GB, it implies different effects of the 
SIAs at the GB on the GB's stability, depending on the GB structure. As seen in Fig. 7(j), the maximal AE,, 
basically increases with E77, linearly. As a result, the loading of the SIA to the GB with a lower SIA stability 


(a large E% 


SIAf 


) would cause the GB to be more easily instable. Figure 7(j) also suggests that A£, exhibits an 
approximately negative relation to the £,,,. Given the meaning of£,,,as the indicator for the intrinsic 
thermodynamic stability of a pure GB and the AZ,, as an indicator for the radiation-stability of the GB, such 


relation may reveal the contradiction of the two types of stabilities. To maintain a relative stability of a GB 


under radiation, a specific intrinsic stability of the GB is required. 

We also found that the profiles in X5(2 1 0)/[0 0 1] (Fig. 7(a)), X3(1 1 2)/[1 1 0] (Fig. 7(e)) and X11(3 3 
2) [1 1 0] (Fig. 7(1)) exhibit a obviously periodical trend similar to that in X5(3 1 Oy[O 0 1] (Fig. 3(a)) 
(labeled as Type 1 in Fig. 7). However, for the other GBs a different trend was observed, where the £,,, firstly 
increases and then basically becomes constant (Fig. 7(b-d) and (f-h) labeled as Type 2). Given the SIA-V 
transition during the GB structural relaxation, the two different types of the GB energy response to the 
accumulated SIAs at the GB could be related to the distinct SIA/V energetic properties at the GB. As shown 
in Fig. 7(j), the GBs of Type 1 were found to have a large absolute difference between CZ; and C7" 
-E®]|) 


(ES 


SIAf 


, while the E77, — E is small in the GBs of Type 2. This could be readily understood, since a 


small value of Dx -Ej implies that the slight variation in the GB energy would be induced after the 
transition of the SIA to V at the GB. 

Based on the structural response of the GB to irradiation, the GBs were further divided into two 
categories of Type I and Type II. The GBs of Type 1 and some of GBs of Type 2 belong to Type I. These 
GBs have similar defects accumulation and structural relaxation behavior to that in X5(3 1 0)/[0 0 1]. The 
loaded SIAs also become components of the GB and new ordered GB structures form (Fig. 7(al—el) and 
(11)). As the concentration of the SIA at the GB increases to a certain level, the GB migrates locally, 
accompanied by the transition of the SIA to V. These GBs were found to have high SIA/V formation energy 
at the GB (E sum) (Fig. 7(k)). The others belonging to Type II exhibit the unique response to the V/SIA 
accumulation at the GBs (Fig. 7(fl-h1)). We frequently observed the formation of the disordered GB 
structures during the loading of the SIAs to the GBs, which could arise from the low thermal stability of the 
GBs as suggested by the small Br in these GBs (Fig. 7(k)). These GBs have a quasi-continuous 
distribution of the E/ after absorbing a certain concentration of SIAs, as shown in supplementary Fig. S13 for 
the case of X9(1 1 4)/[1 1 0]. 


Figure 7(a2—12) further presents the evolution Of Ches; ree2ree3 With C; in the two types of GBs. For the GBs 


of Type I, the similar profiles were observed to that for X5(3 1 0)/[0 0 1], including the alternating variation 


inthe Cre aNd Crez- It indicates the mutual transition of the GB's role of the trapping and annihilation for the 


V during the SIA loading to the GB. However, for the twin boundary of 23(1 1 2)/[1 1 0], we frequently 


andC 


observed the simultaneous increase of the C "T 


T This is due to the very small segregation energy of 


the V at the pure GB (0.18 eV in supplementary Table S8). Consequently, C, increases from zero point. 


Reg 
Therefore, an additional picture about defects accumulation and GB structural relaxation works in the twin 
boundary via the synergistically enhanced trapping and annihilation of the V. For the GBs of Type IL, C, 


egl 


is much larger than Cres and C pegs , indicating that the trapping is a main process for the GB absorbing the V 


eg3 ? 


is close to C, ,; besides, both C, and C 


nearby. It can also be seen that C Reg3 ? Reg2 


Reg? do not change apparently 


Reg3 
with C2, , compared to that in the GBs of Type I. Such variation feature well corresponds to the disordered 


GB structure of Type I. The level of Cpe; may be related to the Ey, in the GBs (supplementary Fig. S12). 


eg3 


The small ES. in the GBs of Type II leads to the small level of C 


neg; » ID agreement with the GB classification 


based on Es, (Fig. 7(k)). 


3.7. OKMC remodeling of defects accumulation and structural relaxation near the GB 

By performing MD simulations of the GB structure predicted by the OKMC simulation and the V-GB 
interaction, we have revealed the new processes of the GB relaxation and its interaction with the V. 
Considering these new atomic processes, we remodeled the defects accumulation and structural relaxation 
near the GB by using the OKMC method. The main features captured into the new OKMC model include 
the disappearing of the SIAs at the GB after a critical Cf’ and the simultaneous GB migration as observed in 
the GBs of Type I. The effect of the formation of the new GB phase on the V nearby is similarly reflected by 
the transition of the annihilation region to the trapping region. As for the GB of Type II, it does not involve 
the SIA-V transition, while the trapping for the V always works (Fig. 7(f2—h2)). In the remodeling of the 
defects behavior, we mainly considered the GB of Type I. Since the GB only directly interacts with the V 


near itself, we thus defined the V concentration near the GB ( C7"). By examining its variation combined 


with the snapshots of the defects microstructure, we explored the evolution of the V-GB interaction 
mechanisms with radiation dose. 

Figure 8 presents the OKMC simulation results of the variation in C, and the corresponding C; ^ with 
radiation dose at several values of T and z (see supplementary Fig. S14 for the results of 7-300 K and r of 
10? dpa:s'). As seen from Fig. 8(a-c), different to the old OKMC model prediction of C% to increase 
continuously with radiation dose (Fig. 2(e-g)), the new OKMC model predicts that C% no longer increases 
after a critical dose at the most considered radiation conditions (Fig. 8(a—c)) except at the condition of T of 
300 K and z of 10? dpa-s! where both the C7; and C7" remain at a low level (blue curves in Fig. 8(c)). This 
is consistent with the MD simulation results and the preset picture in the OKMC model. Noticeable feature 
was also observed on C" ^, At T of 100 K and c of 10° dpa:s'!, it was found that C7"? firstly increases and 
then decreases from 0 to about 0.1dpa (if there is no SIA-V transition process, C; ^ would always decline 
due to the annihilation of the V with the sufficient amount of SIAs collected at the GB) . During this period 
(as illustrated in Fig. 8(a) by P), the SIAs gradually accumulate at the GB (Fig. 8(d) and (e)), and finally a 
region free of Vs but with a high concentration of SIAs at the GB was formed (Fig. 8(e)). The period well 
corresponds to the trapping-annihilation transition stage (from 70 to P/ in Fig. 5(c)). Since 100 K is too low 
to activate the trapping, the annihilation plays a more and more important role in removing the V near the 
GB as more SIAs are stored at the GB. Such increase-decrease trend in C7" ^ was repeatedly observed as 
radiation continues; similar pattern in the microstructure also occurred repeatedly (Fig. 8(d), (e) and (f), (g)). 
Besides, the peak value of C7" becomes large in the following several periods. In addition, there is no such 
V accumulation stage at relative high T of 200 and 300 K (Fig. 8(b) and (c)), which could be attributed to 
the effective absorption and annihilation of the V by the GB due to GB-enhanced diffusion and annihilation 
of the V nearby. Furthermore, after the critical C?» , C^ ^ immediately reaches a peak value. Examining the 
microstructure, we found that these features could arise from the V accumulation in the bulk region near the 


GB (Fig. 8(e) and (g)) and the formation of V clusters (Fig. 8(h-k)). Consequently, after the SIA transition 


and the simultaneous GB migration, more Vs enter the GB interacting region. 


4. Discussion 


4.1. The pattern for the defects accumulation and GB structural relaxation under continuous radiation 

As mentioned in the introduction, the real radiation conditions (T, t, rd, L) are complex and span a large 
parameter range over which radiation damage could accumulate or be healed through distinct mechanisms. 
We summarized the main parameters in irradiation experiments on typical Fe and W systems (supplementary 
Table S1), given the importance of the two systems in the nuclear engineering application. Fig. 9(a’) shows 
the experimental parameters distribution in the T-t space. In the parameter regime of a small L in NCs, a 
high T and a low c (Regime A in Fig. 9(a')), OKMC simulations (Fig. 2, supplementary Figs. S7 and S8) 
have shown that radiation damage could not accumulate in the system. The self-healing property for 
radiation damage in the parameter regime could be accounted for by the basic segregation and annihilation 
processes promoted by the GB, e.g. V-hop and /E [] (of course, IE has also been shown to work at a T [] in 
Regimes B/ and B2). 

In the other regimes of low T and/or high rt) of the authors' interest (Regimes B/, B2 and B3), the 
present study however reveals an alternative pattern or mechanism for the defects accumulation and GB 
structural under continuous radiation, involving the cooperative relaxation of the defects and GB. As 
illustrated in Fig. 9(a-j), the defect-GB interaction evolves with the distribution of the defects or radiation 
dose (Fig. 2) in terms of the change of the energetic and kinetic properties of the V near the GB and at the 
GB, and also the modification of the GB structure, energetics, and the kinetics by the segregated defects at 
the GB (Figs. 3-7, supplementary Figs. S4-S6, S9, S10, S12). Prior to irradiation, a pure GB could serve as 
a sink for the SIA and V due to the well-known energetic and kinetic reasons (supplementary Figs. S5 and 
S6, []), as shown in Fig. 9(a) for the trapping region for the V near the GB. After the SIA escaping the 


annihilation with a bulk V and segregation at the GB, the annihilation region forms around the SIA [] (Fig. 


9(b)); the V near the region could be recombined via previously reported JE at a low energy barrier (process 
(D in Fig. 9(b)) []. Due to the preferential absorption of the SIA over the V, there is a net increase in the SIA 
concentration at the GB (Fig. 9(c)), while the V in the bulk region far away from the annihilation region 
accumulate in the grain interior (Fig. 2(k-n), Fig. 8(d—k)). 

Above a critical SIA concentration, the SIAs at the GB lost their identity accompanying the formation of 
a new GB phase forms (Fig. 9(d)). Near the new GB structure, the V diffusion is further accelerated 
compared to that near a pure GB in Fig. 9(a) (Fig. 6), leading to the trapping of the V by the GB via a low 
energy barrier process @) (Fig. 9(d)). With the relaxation of the defects near the GB and the SIAs at the GB, 
the V in the bulk region far away from the GB meanwhile gets clustered (Fig. 9(d)) or assemble in the 
neighboring hood of the GB. From Fig. 9(a) to (d), the GB energy increases with the SIA accumulation at 
the GB (Fig. 9(1)). During the period, the functional region for the V near the GB evolves from the trapping 
region for the GB (Regimel in Fig. 9(j)) to the annihilation region (Regime3 in Fig. 9(j)) and finally to 
Regime2 with E/ close to zero (Fig. 9(j)). With the SIA continuing to accumulate at the GB, the GB begins 
to recover to its original structure via the transition of the excess SIAs trapped at the GB to the minority of 
the V at the GB and the GB migration (Fig. 9(e)). The GB energy also decreases remarkably, while the 
functional region returns to the trapping region for the V. Consequently, the V on the one side of the GB (e.g. 
V3 in (Fig. 9(e))) that initially locates in the bulk enters the GB’s interaction range and could be trapped by 
the GB. Particularly, the V near the local GB structure containing Vs could also easily segregate to the GB 
(Fig. 9(f)) due to the enhanced binding of the V to the GB and diffusion by the Vs at the GB (Fig. 6). The 
subsequently segregated SIAs could annihilate the V at the GB (process © in Fig. 9(g)). Meanwhile, the 
bulk V-cluster (Vn) located near the GB could also be trapped by the GB due to the GB motion (process ©) 
in Fig. 9(g)). Such cycle (M1 in Fig. 9) repeats over radiation dose (supplementary Table S1) involving the 
collection of the SIAs at the GB, the formation of new GB phase, the SIA-V transition and the GB motion, 


and the interaction of the V/V, with the GB. Depending on the GB structure, the alternative response 


mechanism for the GB cumulative irradiation exists (M2 in Fig. 9). The GB with low defect formation 
energy could also develop to the disordered structure with the defects accumulated within the GB. 

The revealed defects accumulation and GB structural relaxation mechanisms (M/ and M2 in Fig. 9) ina 
model system of Fe would apply to other systems, e.g. NCs Au and W, although the exact transition kinetics 
may depend on the material type and GB characters. The systems of Fe and Au have similar relative defects 
properties, e.g. the higher mobility of the SIA compared that for the V, which induces the common 
phenomenon of the preferential segregation of the SIAs at the GB over the V. Despite the different levels of 
defects migration energy barriers in bulk Fe and W (0.63 versus 1.80 eV for V migration), only the T range 
of a low T is affected. For example, the T of 563 K in irradiation experiments on W [] could be considered as 
a low T for W. These results also motivate investigating the defects accumulation processes and 
interface/boundary relaxation mechanisms in nano-composites [] and oxide dispersion strengthened 


materials []. 


4.2. Implication of the GB relaxation mechanism for designing the radiation resistance of a nano-structured 
metal 

Reducing the grain size in NCs is a well-known option for enhancing the radiation resistance of a 
poly-crystal metal. Our results provide another atomic insight into radiation damage reduction in NC metals 
from the view of the coupling relaxation of the defects and the GB structure, and have implications for 
designing and optimizing radiation-resistant NCs. The present work reveals several compromise approaches 
to tailor radiation performance of NCs. 

First, although the radiation damage near the GB could be healed via the interaction of the defects with 
the GB, the damage at the GB should not be neglected. Generally, it is considered that the defects segregate 
to the GB in NCs and then are eliminated at the GB via complex annihilation mechanisms near the GB []. 
The trapping of the defects by the GB alleviates the defects accumulation in the grain interior, reducing 


possible swelling. Yet, depending on the radiation conditions and the energetic and kinetic properties of the 


defects in the grain interior, near the GB and at the GB, the SIA/V could accumulate within the GB. The 
SIA-V transition mechanism revealed in the present work could prevent the over-accumulation of the 
defects at the GB. Consequently, the GB owes a sustainable capacity for absorbing defects from the grain 
interior. The mechanism is related to the defects behavior within the GB and the GB relaxation mechanism, 
e.g. GB migration. The GB motion leads to competing effects related to the radiation performance. In 
addition to the trapping of the V/SIA by a static GB via the motion of the defects, a moving GB could more 
effectively absorb the V/SIA in the grain interior []. And consequently, the defects accumulation rate in the 
grain interior is lower than that near a static GB. Nevertheless, it is well known that the GB migration would 
lead to the grain growth. The reduction in the trapping efficiency caused by the increase in the grain size 
could cancel the enhanced trapping of the V/SIA by a moving GB. 

Second, there should be compromise between radiation stability of the GB and the intrinsic stability of a 
pure GB. The present study suggests that the GB stability is reduced as radiation defects accumulate at the 
GB. Furthermore, the magnitude of the increase in the GB energy is basically inversely correlated with the 
GB energy of a pure one. 

Third, there is contradiction between the GB absorption rate for defects and the GB's thermal stability. 
To maintain the basic mechanical property of a GB, it is rational to select the GB that has a low V/SIA 
accumulation rate as the one with better radiation performance (Regime A in Fig. 9(a')). To improve the 
trapping efficiency of the defects by the GB, it qualitatively requires the increase in the binding energy of 
the defects with the GB. For example, the twin boundary binds with the V/SIA weakly compared to other 
GBs, and consequently the defects accumulation rate in the grain interior is larger than in other GBs at 400 
K. To obtain better radiation resistance at a high 7, the fraction of the twin boundary in the nano-structured 
system should be decreased. At a low T, e.g. 300 K, the twin boundary behaves similarly to other GBs in 
defects accumulation behavior. Yet, we should not overemphasize the value of the V/SIA-GB binding 


energy. The thermal stability and defect-GB binding strength (binding energy) are two contradictory factors. 


On the one hand, at an elevated temperature, as the binding energy 1s higher than a certain level, the defects 
accumulation in the grain interior is insensitive to the binding energy. On the other hand, a large binding 
energy corresponds to small defect formation energy at the GB. As the defect formation energy is too low, 
the GB exhibits a poor thermal stability. Consequently, as defects accumulate within the GB, the disordered 
structure is easily formed. Such structure may serve as the nucleation point for the radiation-induced 
amorphization []. 

Finally, the present study may also touch the controversy over the equivalence between ion and neutron 
irradiation. As NC Fe is used as structural materials for next-generation fission and fusion reactors, it should 
be able to survive up to 200 dpa of damage at a high operating T of ~550°C []. Yet, the t is as low as 
105—107 dpa-s'. And thus, the relevant parameter regime may locate within Regime A (points A and B in 
Fig. 9(a")). However, in many accelerated irradiation testing experiments (e.g. the use of ion irradiation as a 
surrogate for neutron irradiation) the t is much higher (~107—10° dpa-s’') than that in reality. Besides, some 
experiments are performed at a low T to elucidate the effect of a specific atomic process. And hence, the 
radiation conditions in ion radiation experiments may locate within Regimes 41, A2 and A3 in Fig. 9(a’). We 
have shown that different defect accumulation and GB structural relaxation mechanisms apply to these 
parameter regimes. Therefore, it may be controversial to extrapolate the results about the defect-GB 


interaction and microstructure evolution from high dose rate to low dose rate. 


4.3. The limits in the current predictor-corrector method 

Complementary to experiments, theoretical modeling/simulation plays an important role in 
understanding materials properties degradation and predicting their performance under irradiation. 
Radiation-induced microstructural evolution is an inherent multiscale process ranging from atomistic to 
macroscopic time-space level []. To model this complex process, a multiscale approach with information 
passing between different scales is often required. The defect energetics, kinetics, and controlling 


mechanisms obtained from atomistic modeling are usually used as inputs of the mesoscale methods, e.g. 


KMC. However, the present study suggests that the feedback from mesoscale simulations to atomistic 
calculations is necessary to accurately or precisely predict the microstructure evolution under cumulative 
irradiation. In other words, the close coupling between the two methods in the neighboring scales is needed. 
This motivates the development of the predictor-corrector method in the current work. 

The limits of the present method exist. There is obvious difference in the scale that the two methods of 
the atomic and coarse-grained simulation methods can access, particularly the space scale. It brings about 
the difficulty in dealing with the KMC-passed large clusters in atomistic simulations. This is important, 
since radiation to the material in real fusion/fission or ion radiation conditions produce both point defects 
and complex defect structures, e.g. V clusters and SIA loops. These clusters may have rather different 
behaviors near the GB from the point SIA/V, affecting the structural relaxation mechanism for the GB with 
defects clusters. In addition, the possible coupling between the cascade and the thermally-activated atomistic 
processes is also hard to deal with in the present simulation framework. Finally, besides for the variation in 
the basic energetic and kinetic properties of the defects with radiation dose, the other parameters e.g. the 
sink strength and its evolution with radiation dose relevant to the high level of simulating techniques are also 


needed. 


4.4. The SIA-V transition in the locally pinned GB 

Generally, NCs exhibit poor thermal stability of and coarsen significantly when heated to 0.3 to 0.5 of 
their melting temperature []. To stabilize NCs, nano-precipitates are often introduced at GBs to reduce the 
GB mobility and hence slow the rate of grain growth []. By pinning several atoms near the GB of 2£5(3 1 
Oy[ 0 01], we modeled the drag of the GB by the nano-precipitates and explored whether the observed 
SIA-V transition phenomenon occurs in a stabilized NC. We found that the curves of the GB energy as a 
function of Cf" for a locally pinned GB (supplementary Fig. S15(a)) exhibit similar trend to that for a 
freely-migrating GB (Fig. 3(a)). The GB energy is only slightly higher as some atoms are pinned due to the 


inadequate relaxation of these pinned atoms. We further examined the GB structure and atom displacements 


fields after relaxation (supplementary Fig. Sl5(b-f). As one atom on the GB plane (atom A in 
supplementary Fig. S15(d)) was pinned, the transition occurred with the grain on one side of the GB 
migrating to the neighboring low energy state. As one atom on one side of the GB (atom B in supplementary 
Fig. S15(e)) was also pinned, the grain on the other GB side migrated during the SIA-V transition. 
Nevertheless, as more atoms (atoms A, B and C in supplementary Fig. S15(f)) were pinned, the atoms 
surrounding the pinned atoms also nearly remain immobile, while the GB parts far away from the pinned 
region migrated during the relaxation. The SIA near the pinned atoms did not transit to the V (supplementary 
Fig. S15(c)), but the other SIAs transited to the Vs or V-SIA-V complexes (supplementary Fig. S15(c)). 
Therefore, we conclude that the SIA-V transition would occur within the GB containing the 


nano-precipitates. 


5. Conclusions 

In this study, we developed the predictor-corrector simulation method and preliminarily realized the 
close coupling between the atomistic and coarse-grained modeling techniques. Using the method, we 
investigated the defects accumulation and GB structural relaxation mechanisms in NC Fe under successive 
irradiation. We revealed the radiation-parameter dependent picture for the defect-GB interactions. Within the 
parameter regime of a high temperature and or low dose rate, defects are well annihilated within the GB, 
rendering the GB system exhibit the self-healing property for radiation damage. The performance is 
explained by the GB-promoted annihilation of the V with SIA via the well-known processes of /E and V-hop. 
Nevertheless, within the parameter regime of a high dose rate and/or a low temperature, the new dynamic 
defect-GB interaction process were revealed involving the cooperated relaxation of the defects and GB. As 
radiation proceeds, the SIAs initially accumulate at the GB, and then become structural components of the 
GB, and finally the GB recovers to its background structure. The recovery involves the transition of a critical 


concentration of the SIAs to the minority of the Vs accompanied by the local GB motion. We revealed the 


enhanced location and migration of the V near the GB with a new phase, containing Vs, making the V-type 
defects instable near the GB. By calculating the cumulative distribution function of the vacancy formation 
energy and the grouping of the sites into several regimes with the same physical meaning, we meanwhile 
revealed the alternating evolution of the GB structure and role of trapping and annihilation of Vs. We also 
found that some GBs with low thermal stability develop into disordered structure during irradiation. The 
results reveal the evolution of the distribution of the defects near the GB, the energetics and kinetics of 


defects, GB energy, structure and mobility with radiation dose. 
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Fig. 1. (a) Geometric model for OKMC modeling radiation damage accumulation in NC Fe. The grain size 
is L. The GB locates at the center of the cubic; planes ABCD, EFGH, IJKL are GBs. (b) Schematic 
illustration of the combination of the MS, MD and OKMC methods simulating defects structure evolution 
under accumulative radiation. The red sphere represents a SIA and the green cubic is for a V. Here, /—9 
represent the fundamental atomic processes incorporated into the OKMC model, including SIA/V diffusion 
and clustering in the bulk (//2), SIA-V annihilation in the bulk (3), IE-induced annihilation of the SIA-V 
near the GB (4), SIA/V segregation from the bulk to the GB (5/6), SIA/V migration and clustering at the GB 


(7/8), and SIA-V annihilation within the GB (9). Csr, y denotes the SIA/V concentration in the simulation 


model. Egg means the GB energy. Actually, the curve of Cs; y- dose illustrates the variation of the SIA/V 
distribution and concentration with the radiation dose. Similarly, the curve of Egg~ Csj4 y denotes the Csr,y 
dependent GBA energy and SIA/V-GB interactions. (c) Illustration of the interaction between two defects in 
the bulk region. Here, the large red sphere represents a SIA-cluster, and the large cubic is for a V-cluster. (d) 
Illusion of the SIA,/V, segregation near the GB. (e) Illustration of the SIA,/V, behavior within the GB. (f) 
shows the energy barrier profiles for the SIA,/V, diffusion and their annihilation near the GB. (g) and (h) 


Illustration of the two styles of producing a certain concentration of the SIA,/V, within the GB. 
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C > Fig. 2. (a-d) The bulk SIA/V concentration (C7/,) as a function of the radiation dose at different 


. temperatures of 100, 200, 300 and 600 K. (e-h) The SIA/V concentration at the GB ( C$, ) as a function of 


the radiation dose at the above four temperatures. The legend for (a-h) is given in (i). In (e-h), the horizontal 


dash line separates the two possible regimes where the GB responses to irradiation via different mechanisms. 


a ; Symbol A in (a) and A/ in (e) respectively mark defects accumulation in the bulk and at the GB under a 
t certain radiation dose. (j) Illustration of the SIA/V concentration level in the bulk and at the GB after about 1 


dpa in the regime of temperature (T) and grain size (L). Here, Hg and Lg respectively represent the high and 


low C" , while Hg and Lg represent the corresponding high and low C, (here, the high level of defects 
concentration means the concentration increases with radiation dose). (k-n) The bulk V concentration as a 
function of the distance of the V to the GB at different temperatures of 100, 200, 300 and 600 K. The legend 
for (k-n) is given in (o). Note that, in (a-h) and (k-n) L is 10 nm. As calculating C/" , the Vs near the GB 
were also considered as bulk defects. In (k-n), the curves were smoothed by using a moving average filter 


with a span of 5. 


202001.00046v1 


chinaXiv 


«— — — — Cycle — ———— «— — — — Cycle1— > «— — Cycle? —_—__> 


1.4-(3)Z5( ) K, Style2 


3GB* 
'SIA2 6 


12 
GB -2 
C i (nm ^) 
Style2:4x4-300 K (or Stylel: 100, 200 and 300 K) 
(e) G2, 7, 7. S CER 
Co On 9, 09, 9 e, 0 


JGB* Jj 
2C 4 c 


B.C 
' e 


[3 1 0] (À) 


rate SB 
Ld 


0 5 10 15 20 25 30 35 5 10 15 20 25 30 35 5 10 15 20 25 30 35 


[1 3 0] (À) 
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1] (CS4) in the manner of Style2 (Fig. 1(h)). Error bars here indicate the variation range of the calculated Ecg 
from multiple independent MD simulations, while the pink dash line indicates the variation trend of the 
profile. Eggo is for the energy of a pure GB without SIAs loaded, while A£., is the variation in Egg relative to 
Eggo. (b-j) The relaxed structures of the GB loaded with several different concentrations of SIAs. The 
horizontal dash black line indicates the initial GB position prior to the structural relaxation of the GB with 
SIAs. In (f-3), the horizontal pink line indicates the GB position after relaxation. In (b-j), the atoms are 
colored with their potential energies (E), as shown by the colorbar on the right hand. In (a), Co is zero, while 


symbols Cp, Ce, Ca, Ce, Ch Ce, Cn, Ci, C; respectively denote the Cj; in (b-j). In (e), the inset figure shows 


the distribution of the atom separation distance in the two GB systems with SIA concentrations of 0.0 and 
about 4.6 nm”. The calculated regions are limited to the atom layers near the GB with the width of 6 A. For 


the pure GB, such region consists of three atom layers on each side of the GB. 
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| . in the manner of Style2. (a) shows the atom motion process during the transition of the GB system by 
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and green lines respectively indicate the initial and final GB positions. The black circle and green polygon 
indicate the GB structural unit before and after transition, respectively. (b) shows the corresponding atoms in 
initial state GBs1(4, B, C, D, E, F) to that in final state GBs2 (41, Bl, CI, DI, El, FI). Note that 
AIIBI/CIIDIIEI/FI mark a specific column of atoms aligned along [0 0 1] in GBs2. The atoms in one 
column are further indicated by /, 2, 3, 4. In GBs1, the corresponding column of atoms A/B/C/D/E/F in 
GBsl are not well along [0 0 1] due to the distortion induced by the loaded SIAs. (c) shows the transition of 
the atoms G in GBs1 to G/ in GBs2. Note that the gray sphere denotes the atom that actually does not exist, 


and hence the gray square represents the vacancy formed in atom column G/. (d) and (e) present the similar 


transition of the GB as the SIA is put in to the GB in the manner of Style] at 0 K. (f) Illustration of the 
SIA-V transition process during the GB relaxation loaded with excess SIAs. /R1, IR2 and /R3 denote the 
atom columns in the initial state, while FR/, FR2 and FR3 represent the atom columns in the final state after 
relaxation. The red spheres denote the interstitials trapped at the initial GB state, while the light green cubic 
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Fig. 5. (a) Cumulative distribution function C( E/ ) in several GB systems loaded with different 
concentration levels of the SIAs. Symbols Reg/, Reg2 and Reg3 represent three regimes with different site 
stabilities; in Regl, E/ < E/ < E" - E,, in Reg2, - E/ « E; « E; , and in Reg3, E/ <- E; . Here, E, and E/' were 
assigned to be 0.3 and 0.2 eV, respectively. (1) shows the average maximal displacement of the atoms after 
the structural relaxation at 0 K of the GB with a vacancy at a specific site in the three regime as a function of 


Cy» represented by dis MaXpegireg2reg3: (J) shows the site density of the three regimes as a function of Cy; 
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( Crceiree2.reg3 ), defined as the number of sites in a specific E/ regime divided by the GB area. Note that, the 
curves in (j) were smoothed by using a moving average filter with the span for the moving average of 5. 


Co is zero, while symbols Cp, Ce, Ca, Ce, Cj, Ce, Ch, Ci, C; respectively denote the Cy; in Fig. 3(b-j). 
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Fig. 6. (a-e) show the energetic and kinetic properties of the V near several typical GB systems with 


different concentrations of SIAs. In (a-d), the sites with different stability are shown by different colored 


v:202001.00046v1 


x symbols. The horizontal dash line indicates the GB position. (e) The energy profiles for the V diffusion 
c along several typical paths near the GB. The two horizontal red and blue dash lines mark zero point and bulk 
© El, respectively. The corresponding paths are illustrated in (a-d). (f-h) Interaction of the V with GB X5(3 1 
0)/[0 0 1] that has distinct structures revealed by MD simulations at 300 K. (f) The GB is pure without SIA 
loading, (g) the GB has a meta-stable phase, and (h) the GB consists of Vs. A vacancy is denoted by a green 
square. The horizontal dash blue line indicates the GB position. The atoms positions before and after 50 ps 
at 300 K were connected to show the dynamic behavior of the V. In the three systems, the Vs were created 
initially at the same sites. The minimal distance of the V from the pure GB is about 5 À. The V 


concentration near the GB is about 7.6 nm". 
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Fig. 7. (aci) GB energy profile as a function of the loaded C3; in the manner of Style2 (Fig. 1(h)). Error bars 


0 510152025 0 5 10 15 20 25 0 5 1015 2023 5 10 15 2023 5 10 15 20 25 
GB -2 m^ GB £. 
C^. (nm Eg, m) C=" (nm 
SIA ( ) 0204 0608101214 SIA ( ) 
ety (TOT p SN c CSRS se eee 1.2 23(112)/[110] s4 
EET TETTEI S 5 SCAND QU fatter ROMS OR 
a 92825 o 89 g BP o. PIC DRE ae 9089 90007 aot Oe CIE = 1.0 
TELLE I pet a t e ER 
d a AGAAT; A 7$ 08 
eS Oa z $- 
= i 06 : 
© EF : 
«t E 
© z 0.2 buts. 
« 
e s 
" .0 0.5 1.0 1.5 2.0 2.5 3.0 
TT EC? (eV) 


here indicate the variation range of the calculated Egg from multiple independent MD simulations. Here, 


AE, 1S E,, referred to the GB energy of a pure GB (Ess). (al-il) The GB structures loaded with a 


naXiv 


c concentration of about 25.0 nm? SIAs. The atoms are colored with their potential energies by using the 


C 


identical scheme to that in Fig. 3. In (a-1i), axis Z is along [2 1 0], [5 1 0], [3 2 0], [1 1 1], [1 1 2], [1 1 4], [2 
2 1], [1 1 3], and [3 3 2], respectively. In (a—c), axis Y is along [1 2 0], [1 5 0], [2 3 0], respectively. In 
(d-i), Y is along [1 1 0]. (a2-i2) The site density of the three regimes Reg/, Reg2 and Reg3 as a function 
of Cg. Symbols Reg/, Reg2 and Reg3 have the same meaning as that in Fig. 5. Here, E, and E/ were also 
assigned to be 0.3 and 0.2 eV, respectively. (J) shows the relation of the maximal AZ,,, to the formation 
energy of a single SIA at a pure GB (fas bottom axis) and £,,, (top axis). (k) shows the absolute difference 
between the SIA and V formation energy at the GB (|z¢", - £67 |) (left axis) and the minimal SIA/V formation 


SIAf 


energy at the GB (right axis) as a function of the GB type. 
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Fig. 8. (a-c) Variation inc% (the left axis) and the corresponding V concentration near the GB ( C7") (the 
right axis) with radiation dose at several values of T and r. According to Fig. 3(a), the critical SIA 
concentration at the GB was set 4.0 nm” for inducing the GB motion. Each time, the GB migrates along a 
random direction jumping over a distance of 2.5 À. The grain size is 10 nm. The region near the GB was 
defined as the atom layers across the GB with a width of W. W was assigned to be 10 À in most cases. At T 


of 200 K and z of 10? dpa: s'', W was given of 16 A. At T of 300 K and t of 10? dpa: s'!, W was given of 30 


A. Ata T of 300 K and t of 10? dpa: s'', W was given of 16 A. At such relatively high T and /or low q, there 
is no defects survived near the region with a certain distance from the GB. Here, C7" was defined as the 
number of Vs in the region near the GB divided by the volume of the region. The volume was calculated as 
W times the GB area. (d-k) shows typical snapshots for the microstructure evolution at an elevated T and z. 
In this visualization scheme, a V-cluster is displayed by a red square, while a SIA-cluster is shown by a red 


filled circle. The radius of the square/circle is proportional to the real size of the V/SIA-cluster. 
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Fig. 9. (a') Illustration of the T-t parameter regime where the data points are obtained from references []. (a—j) 


Schematic illustration of two kinds of the GB response patterns to cumulative irradiation. A SIA is denoted 


by a red sphere and a V is represented by a green cubic. The light blue sphere indicates the atoms at the GB, 


illustrating the GB structure and acting reference atoms for the SIA/V site within the GB. The vertical dash 


line in (a) indicates the GB position. The pink filled ellipse in (b), (c) and (g) represents the annihilation 


region formed around a SIA trapped at the GB. The deep blue filled ellipse in (d) represents the extended 


enhanced diffusion region for V near the meta-stable GB phase, while the large blue sphere represents the 


SIA that loses its identity at the GB. Symbols W0, W1 and W2 respectively represents the width of the pure 
GB, the GB with SIAs, and the meta-stable GB for trapping/annihilating the V. Symbols Eo, E1, E>, E; and 
Eo’ respectively denotes the energy for the pristine GB, the GB with small concentration of SIAs, the GB 
with high concentration of SIAs, the meta-stable GB and the GB containing a certain concentration of Vs. 
symbol dp in (e-g) denotes the distance for the GB migration after the transition from (d) to (e). (i) illustrates 
the energy and structural variation during a response period of the GB to cumulative irradiation. Symbol © 
marks the JE-induced annihilation of the SIA with the V near the GB, (2) marks the trapping of the V by the 
meta-stable GB, (3) marks the trapping of the V by the GB containing Vs, (4) marks the trapping of the V by 
a migrating GB, © marks the annihilation of the V at the GB with the SIA nearby, and © denotes the 


clustering of the V near the GB. 
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1. Supplementary figures 
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: ‘ Fig. S1. The number of bulk vacancies (N;"" ) after the irradiation to the symmetric tilt grain boundary (GB) 
= X5(3 1 0)/[0 0 1] by many consecutive 3 keV-PKAs (primary knock-on atom, PKA) at a time interval of 50 
= ps in a-iron (Fe). Results from Ref. [S1] are also shown for comparison. Molecular dynamics (MD) was 
| e : performed to study cascade-induced damage near the GB at 1000 K using the velocity-Verlet method for the 
© numerical integration, and the simulations ran for about 5000 ps. An atom at one side of the GB was 
selected as the PKA. The atom was given 3 keV of kinetic energy with its velocity directed perpendicularly 
toward the GB. The atoms in the outermost layers of the moving region, with a thickness of two times the 
lattice constant, were coupled with a velocity-rescaling thermostat to absorb the cascade energy and 
maintain the system temperature at 1000 K. The embedded-atom-method (EAM) interatomic potential 


developed by Mendelev et al was used to model the interatomic interaction [S2], and the short-range form 


of the potential was splined to reproduce the high-energy empirical potential of Ziegler et al [S3]. 
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Fig. S2. Snapshots for the accumulation of vacancies nearby the GB and their absorption by the migrating 


naXiv 


c GB. (a-g) The growth of accumulated vacancies in the bulk region during the successive radiation to the GB 


C 


with the PKA randomly selected near the GB. (h—l) The dissolution of accumulated clusters and absorption 
by the migrating GB. The arrow indicates the vacancy position. V denotes the mono-vacancy. The dashed 
circle sketches out the vacancy clusters. The axes X, Y and Z are along [310], [ 1 30] and [001], respectively. 


The atoms within the GB are colored with their potential energies, while the defects near the GB are red. 
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Fig. S3. Schematic illustration of the factors determining the radiation damage accumulation in the bulk/GB 
region of a nano-crystal. These factors include radiation conditions of temperature (7), radiation dose (rd), 
radiation dose rate (t), and the grain size (L). The atomic processes in the NC system include the processes 
in the bulk represented by Pyu(V, SIA), the processes near the GB denoted by Paeargg( V, SIA), the 
processes within the GB represented by Patap(V, SIA), the process associated with GB structural relaxation 
denoted by Pas(GB structure). The GB region is separated from the bulk region by two horizontal dash red 


lines. The GB plane is indicated by a horizontal red line. 
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= Fig. S4. The relaxed structures for the investigated thirteen tilt symmetric GBs. The atoms are colored with 
~ their potential energies. (a-g) are for the GBs with tilt axis of [0 0 1], while (h-m) show the GBs with the 
C 

c tile axis of [1 1 0]. In (a—m), axis Z is along [2 1 0], [3 1 0], [5 1 0], [32 0], [A3 0], [7 6 0], [13 1 0], [1 1 1], 
io, 


[1 12], [1 1 4], [2 2 1], [1 1 3], and [3 3 2], respectively. And axis Yis along [1 2 0], [1 3 0], [1 5 O], 
[2 3 OL[3 4 OL[6 7 0], [1 13 0], respectively in (a-g). In (h-m), Y is along [1 1 0]. In each figure, 


the tilt angle and the GB energy are marked. The horizontal dash line indicates the GB position. 
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Fig. S5. (a) SIA formation energy ( E$} ) profile near the GB of $5(3 1 0)/[0 0 1] as a function of the 


Distance of V to ih (À) [3 l 0] 


distance of the initial SIA center to the GB. (b) SIA configuration within a GB structural unit, as labeled by 
symbols 77 and 72. A SIA is denoted by a red sphere. The EZ, is marked near the SIA. (c) V formation energy 
(EL) profile near the GB as a function of the V distance to the GB. (d) V configuration within a GB 
structural unit. A V is represented by a green square. V1, V2, V3 and V4 denote the four V occupation sites. 
Ef is marked near the V. In (a-d), the dash black line indicates the GB position. In (b) and (d), the polygon 


indicates the GB structural unit. (e) Fundamental parameters charactering the SIA/V behavior near the GB, 


including SIA formation energy at the GB ( EC^ 


sar). SIA segregation energy (£77), SIA emission energy 


seg 


barrier from a pure GB ( EF), V formation energy at the GB ( Ep Os V segregation energy ( E". ), V emission 


seg 


energy barrier from a pure GB ( E^), the migration energy barrier for the SIA along the tilt axis ( E), ), the 


SIA migration energy barrier in the direction nonparallel to the tilt axis ( E), the migration energy barrier 


im 


for the V along the tilt axis ( E, ), the V migration energy barrier in the direction nonparallel to the tilt axis 


vm 


( E! ), binding energy of the V-V ( E; ' ) and SIA-SIA ( £7" ^") at the GB. The values for the corresponding 


vm 


parameter in the bulk are also given for comparison. 
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Fig. S6. (aj) The V/SIA configurations at several GBs. A V is shown by a green square; a SIA is 
represented by red spheres, the number of which depends on the SIA configuration at the GB. For the V, the 
metastable site is also shown with V formation energy above the minimal E/ less than 0.2 eV. The E/ is 


marked on a site. E/, is marked near the SIA configuration. The vertical blue line indicates the GB position. 


SIA 


In (1) and (k), the inset shows the SIA configuration along another view. Axis Z is normal to the GB plane, 


while axis Y is parallel to the GB plane. Axes Y and Z are identical to that in Fig. S4. 
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Fig. S7. (a-l) The bulk SIA/V concentration (C5,,) as a function of the radiation dose at different 
temperatures of 100, 200, 300 and 600 K and grain sizes of 50, 100 and 200 nm. (al-11) The SIA/V 
concentration at the GB (C5, ) as a function of the radiation dose at the above four temperatures and grain 


sizes. The legend for (a-l) is given in (m). The legend for (a1-11) is given in (m1). 
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Fig. S8. (a-l) The bulk V concentration as a function of the distance of the V to the GB at different 
temperatures of 100, 200, 300 and 600 K and dose rates of10?, 10? and 10? dpa-s''. The legend for (a-l) is 


given in (a). 
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Fig. S9. The relaxed GB structures of $:5(3 1 0)/[0 0 1] with a certain concentration of SIAs at 300 K for 50 
ps. In (a) and (b), the horizontal black line indicates the GB position. The polygon indicates the GB 


structural unit. The atom is represented by a sphere colored by the potential energy of the atom. 
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Fig. S10. (a) The trajectories of a SIA near all the eleven GBs at 300 K drawn from object kinetic Monte 
Carlo (OKMC) simulations of the SIA behavior near the GB. Initially, the SIA is about 6 nm away from the 
GBs. The red curve is for &3(1 1 2)/[1 1 0]. The horizontal blue line indicates the GB position. (b) The 
trajectories of a V near all the GBs at 300 K. (c) and (d) shows the trajectories at 400 K for a SIA and V, 


respectively. 
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Fig. S11. The GB migration distance within 50 ps at 300 K in several distinct MD simulations, which is in 


unit of the distance that the GB moves from one stable state to another nearest one (d0). 
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Fig. S12. The relation of the Creg3 to the formation energy of a single SIA at a pure GB (Egy ). 


202001.00046v1 


chinaXiv 


29(114)/[110] 


| rà 

—n— C =0.00 
GB, 

—o— CC 20.24 


GB 
—^— CH =0.48 


GB 
—7— CSF -9.72 


—0— CH" =0.96 


-1.2 -0.8 -0.4 0.0 0.4 0.8 1.2 1.6 2.0 


Fig. S13. Cumulative distribution function C( E/ ) in the GB X9(1 1 4y[1 1 0] loaded with different 


concentration levels of the SIAs. 
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) Fig. S14. (a) Variation inC, (the left axis) and the corresponding V concentration near the GB ( C7/"^^) (the 
. right axis) with radiation dose at a temperature (T) of 300 K and dose rate (7) of 10? dpa-s”’. (b-e) Snapshots 
for microstructure evolution at above 7, t condition. In this visualization scheme, a V-cluster is displayed by 


‘a red square, while a SIA-clusters is shown by a red filled circle. The radius of the square/circle is 


: proportional to the real size of the V/SIA-cluster. 
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Fig. S15. (a) GB energy (E, ) as a function of the SIA concentration at the GB of X5(3 1 0)/[0 0 1] (CS). (b) 
and (c) The initial structure of the GB with a certain concentration of SIAs and the GB structure after 
relaxation at 300 K, respectively. Three atoms near the GB were pinned during the relaxation. The atoms are 
colored with their potential energies, using the same visualization scheme as that in Fig. 4(c—h). (d-f) The 
GB relaxation process. A, B and C mark the pinned atoms during the structural relaxation. In (f), the atoms 


that nearly remain immobile are covered by the light green circle. In (d—f), the black and blue horizontal line 


indicates the GB position before and after relaxation, respectively. 


2. Supplementary tables 


Table S1. Irradiation conditions in several systems. 


Reference 
[S4] 
[S4] 
[S5] 
[S6] 
[S7] 
[S8] 
[S8] 
[S9] 
[S10] 
[S11] 
[S11] 
[S12] 
[S13] 
[S13] 
[S14] 
[S14] 
[S15] 
[S15] 
[S16] 
[S17] 
[S18] 
[S19] 
[S20] 
[S21] 
[S22] 
[S22] 
[S23] 


[S24] 


System 
Pd 
ZrO, 
Pt 
Ag 
Ni 
Au 
Au 
bcc-Fe 
bcc-Fe 
bec-Fe 
bcc-Fe 
bec-Fe 
bcc-Fe 
bcc-Fe 
fcc-Fe 
fcc-Fe 
fcc-Fe 


fcc-Fe 


Radiation-particle 
4MeV-Kr^* 
4MeV-Kr^* 

2.8MeV-Au* 
1MeV-Kr^* 
1MeV-Kr^* 
60MeV-"C 
60MeV-"C 
100keV-He?* 
10keV-He^* 
lO0keV-He^* 
1O0keV-He^* 
440ke V-He?* 
3.5MeV-Fe* 
3.5MeV-Fe* 
1MeV-Kr^ 
1MeV-Kr^ 
3.5MeV-Fe* 
3.5MeV-Fe* 
4.5MeV-Cu^* 
1MeV-Kr^* 
200keV-He?* 
200keV-He^ 


SkeV-He^* 


3MeV-Si/Cu /W** 


3MeV-Cu^* 
3MeV-Cu^* 
2keV-He? 


1MeV-Kr^ 


T(K) 
298 
298 


573 


298 
298/1050 
298/1050 

1223 


298/1073 


rd (dpa) 
170-210 


3-8 


unknown 

1.03-5.72 
0.2-4 
0.2-4 

unknown 


2 


t (dpa:s )) 
3x10"! 
9x 107 
8x10* 

1.25x10? 
3x10? 

unknown 
unknown 
8x10* 
1x10? 
1x10? 
1x10? 
4x10° 
1x10? 
1x10? 
2.5x107 
2.5x107 
3x10? 
3x10? 
unknown 
0.43-1.98x10° 
2.8x107 
unknown 


unknown 


6x107/9x107/3x107 


1x10* 
1x10? 
unknown 


2x10? 


L (nm) 
10-300 
10-300 
20-100 
100 
200 
23 
23 
49-96 
100 
50 
50 
12.6 
320 
2000 
45 
45 
100 
35000 
280 
unknown 
15 
25-3900 
40-300 
60-400 
50.6-800 
50.6-800 


60-100,100-500 


35,85,80-400,1-3000 


Table S2. The migration energy barriers for small V,/SIA, in the bulk and near the GB. 


Value 


0.67 eV 


0.62 eV 


0.35 eV 


0.48 eV 


2.8eV 


0.34 eV 


0.43 eV 


0.43 eV 


0.43 eV 


2.8eV 


0.0 eV 


0.34 eV 


0.31 eV 


0.18 eV 


0.24 eV 


1.4eV 


Term 


Em_V_bulk 


Em V2 bulk 


Em V3 bulk 


Em V4 bulk 


Em Vn bulk 


Em SIA bulk 


Em SIA2 bulk 


Em SIA3 bulk 


Em SIAnl bulk 


Em SIAn2 bulk 


Em SIAn nearGB 


Em V nearGB 


Em V2 nearGB 


Em V3 nearGB 


Em V4 nearGB 


Em Vn nearGB 


Comments 


Migration energy barrier for a single V in Fe bulk 


Migration energy barrier for V; in Fe bulk 


Migration energy barrier for V3 in Fe bulk 


Migration energy barrier for V4 in Fe bulk 


Migration energy barrier for V, (n>4) in Fe bulk 


Migration energy barrier for a single SIA in Fe bulk 


Migration energy barrier for SIA; in Fe bulk 


Migration energy barrier for SIA; in Fe bulk 


Migration energy barrier for SIA,; (n123, n1«—10) in Fe bulk 


Migration energy barrier for SIA» (n2>10) in Fe bulk 


Migration energy barrier for SIA, (n<=500) near Fe GB 


Migration energy barrier for a single V near Fe GB 


Migration energy barrier for V; near Fe GB 


Migration energy barrier for V3 near Fe GB 


Migration energy barrier for V4 near Fe GB 


Migration energy barrier for V, (n>4) near Fe GB 


Table S3. Some parameters describing the SIA migration along the GB, including migration energy barrier 


for the SIA along the tilt axis (£',), the time for SIA migration along the tilt axis at 300 K (t,), the SIA 


im um 


migration energy barrier in the direction nonparallel to the tilt axis (£/,), the time for the SIA migration in 


im 


the corresponding direction at 300 K (zt). 


im 


GB E; (eV) tn (s) E} (eV) i}, (s) 
bulk 0.33 3.5x107 0.33 35x10" 
£5(2 1 0)/[0 0 1] 0.12 1.0x10 7? 1.52 3.4x10P 
X5(3 10)/[0 0 1] 0.44 25x10" 1.176 57x10" 
X13(5 1 0)/[0 0 1] 0.49 1.7x10^ 3.06 2.6x10? 
x13(3 2 0)/[0 0 1] 0.04 4.7x10 3.05 LTXIQ" 
x25(4 3 0)/[0 0 1] 0.52 5.4x107 1.80 L7*10" 
X85(7 6 0)/[0 0 1] 0.77 8.6 1.7 3.6x10!6 
X85(13 1 0)/[0 0 1] 0.34 52x10" 2.0 4.0x10?! 
x3(11 1)/[1 10] 0.69 3.9x10" 0.50 25x10~ 
©3(1 1 2)/[1 10] 0.55 17x10? 0.07 1.5x10 
x9(114)[1 1 0] 0.45 3.6x10? 1.04 3.0x10? 
xo( 2 1)/[1 10] 0.53 8.0x 107 0.82 6.0x10! 
Xl1(113)[I 1 0] 0.66 1.2x10" 1.22 3.1x10? 


£11(3 3 2)/[1 1 0] 0.87 4.1x10? 0.99 4.3x10* 


Table S4. Some parameters describing the V migration along the GB, including migration energy barrier for 


the V along the tilt axis ( E, ), the time for V migration along the tilt axis at 300 K (t!,), the V migration 


vm vm 


energy barrier in the direction nonparallel to the tilt axis (£",), the time for the V migration in the 


vm 


corresponding direction at 300 K ( ;?. ). 


vm 


GB E! (eV) t, (S) E} (eV) il (s) 
bulk 0.63 3.8x10? 0.63 3.8x10° 
£5(2 1 0)/[0 0 1] 0.36 1.1x10% 1.07 9.5x10° 

X5(3 1 0y/[0 0 1] 0.55 1.7x107 0.82 6.0x10! 
X13(5 1 0)/[0 0 1] 0.68 2.7x10" 1.46 3.4x10" 
X13(3 2 0)/[0 0 1] 0.37 1.6x10° 1.61 1.1x10"° 
x25(4 3 0)/[0 0 1] 0.72 1.2x10! 1.28 3.2x10? 
X85(7 6 0)/[0 0 1] 0.78 1.3x10! 14 3.3x10!! 
X85(13 1 Oy[0 0 1] 0.89 9.0x10? 1.73 1.2x10" 
x3(11 1)/[1 10] 0.95 9.1x10? 0.95 9.1x10? 
©3(1 1 2)/[1 10] 0.60 1.2x10? 0.38 2.4x10° 
r9(1 1 4)/[1 1 0] 0.68 27510" 1.75 2510" 
xo( 2 1)/[1 10] 135 1.0x10? 1.73 12x10" 
L11(1 1 3)/[1 1 0] 0.79 1.9x10! 2.18 4.2x10”% 


£11(3 3 2)/[1 1 0] 0.66 12x10! 0.56 2.6x10? 


Table S5. The binding energy for the V/SIA with small V,/SIA, in Fe bulk. 


Value Term Comments 

0.30 eV Eb V V bulk Binding energy for the V with a V in Fe bulk 

0.37 eV Eb V V2 bulk Binding energy for the V with a V2 in Fe bulk 

0.62 eV Eb V V3 bulk Binding energy for the V with a V3 in Fe bulk 

0.62 eV Eb V Vn bulk Binding energy for the V with a Vn (n>3) in Fe bulk 
0.80 eV Eb SIA SIA bulk Binding energy for the SIA with a SIA in Fe bulk 

0.92 eV Eb SIA SIA2 bulk Binding energy for the SIA with a SIA2 in Fe bulk 

1.64 eV Eb SIA SIA3 bulk Binding energy for the SIA with a SIA3 in Fe bulk 

1.64 eV Eb SIA SIAn bulk Binding energy for the SIA with a SIAn (n>3) in Fe bulk 


Table S6. Binding energy of the V-V ( E; " ) and SIA-SIA ( £7" ^") at the GB. Note that, before calculating 
the binding energy, the stable sites for V/V; and SIA/SIA; at the GB were found by relaxing the 


corresponding GB system at 300 K. 


GB E,” (eV) Ey eV) 
bulk 0.30 0.80 
Xs(210)y[00 1] 0.44 -0.27 
Xs(310)[0 0 1] 0.55 0.25 
X13(5 1 0y[0 0 1] 0.42 0.38 
X13(3 2 0y[0 0 1] 0.52 -0.06 
X25(4 3 0y[0 0 1] 0.63 0.42 
X85(7 6 0y[0 0 1] 0.57 0.40 
X85(13 1 0)[0 0 1] 0.53 0.19 
X3(11 1)/[1 1 0] -1.09 0.02 
£3(1 1 2)/[1 1 0] 0.1 0.48 
X9(114y[110] 0.15 -0.23 
X9(2 2 1)/[1 1 0] -0.25 0.09 
Z11(1 13V[1 10] -0.08 -0.04 


X11 3 2)/[1 1 0] 0.44 0.07 


Table S7. Some parameters describing the interaction of a SIA with a GB, including the minimal SIA 


sur) near the GB, the maximal SIA segregation energy (E; ), the maximal SIA 


formation energy (E 
emission energy barrier from the GB ( E), the time for SIA emission from the GB at 300 K (t), the 


minimal di-SIA formation energy near the GB (E a 


). Note that, the GB with a SIA or di-SIA was relaxed 


at 300 K to find the stable site for the SIA/SIA; location at the GB. 


GB Eg, (eV) Ex (eV) E; (eV) C Ess (eV) 
bulk 3.53 0.00 0.33 3.5x107 
X5(210)/[0 0 1] 0.64 2.89 3.22 1.3x10” 1.55 
X5(3 1 0)/[0 0 1] 0.85 2.68 3.00 2.9x10?* 1.45 
X13(5 1 0)/[0 0 1] 0.71 2.82 3.15 8.8x10?? 1.04 
x13(3 2 0)/[0 0 1] 0.66 2.87 3.20 6.1x10*! 1.38 
x25(4 3 0)/[0 0 1] 0.73 2.80 3.13 3.5x10?? 1.04 
X85(7 6 0)/[0 O 1] 0.59 2.94 3.27 5.8x10^! 0.78 
X85(13 1 0)/[0 0 1] 0.56 2.97 3.30 25x10" 0.93 
x3(11 1)/[1 10] 0.25 3.28 3.61 3.9x10"8 0.48 
©3(1 1 2)/[1 1 0] 2.82 0.71 1.04 2.8x10? 5.16 
x9(114)[1 10] 0.52 3.01 3.34 1.3x10** 127 
xo( 2 1)/[1 10] 0.40 3.13 3.50 1.3x10*6 0.71 
L11(1 1 3)/[1 1 0] 0.46 3.07 3.40 1.410" 0.96 


¥11(3 3 2)/[1 1 0] 0.93 2.60 2.93 1.5x10°7 1.79 


Table S8. Some parameters describing the interaction of a V with a GB, including the minimal vacancy 


formation energy (ES ) near the GB, the maximal V segregation energy (Eke) the maximal V emission 


energy barrier from the GB ( E/^), the time for V emission from the GB at 300 K (t,;), the minimal di-V 


formation energy near the GB ( Ej. ). Note that, the GB with a V or di-V was relaxed at 300 K to find the 


stable site for the V/V» location at the GB. 


GB Ey (eV) E; (eV) E (eV) ty; (s) Ej, (eV) 
bulk 1.71 0.00 0.63 3.8x10? 
X5(2 1 0y[0 0 1] 1:22 0.49 1.12 5.6x10* 2.00 
X5(3 1 0y[0 0 1] 1.20 0.51 1.11 3.8x10* 1.85 
x13(5 1 0)/[0 0 1] 0.85 0.86 1.49 9.3x10" 1.28 
x13(3 2 0)/[0 0 1] 0.80 0.91 1.54 8.2x10P 1.08 
x25(4 3 0)/[0 0 1] 0.78 0.93 1.56 1.7x10'4 0.93 
X85(7 6 0)/[0 0 1] 0.66 1.05 1.34 3.3x10'° 0.75 
X85(13 1 Oy[O 0 1] 0.71 1.00 1.63 2.6x10? 0.89 
x3(11 1)/[1 1 0] 0.35 1.37 1.99 3.2x10?! 1.79 
x3(112)y[l 1 0] 1.53 0.18 0.81 3.8x10! 2.96 
X9(11 4)/[1 1 0] 0.70 1.01 1.62 1.9x10"° 1.25 
x9(22 1)/[1 1 0] 0.34 1.37 2.00 3.3x10?! 0.93 
L11(1 13)/[1 1 0] 0.55 1.16 1.79 1.3x10'8 1.18 
Y11(3 3 2)/[1 1 0] 1.14 0.57 1.20 1.6x10* 1.84 


3. Supplementary text 


Text S1. Derivation for the linear increase of the GB energy with the defects concentration at the GB. 


The meaning of the relevant symbols is defined as follows. S,,: the area of a GB, N : the number of 
atoms in a GB system, £,,,: cohesive energy per atom in the bulk, n: the number V/SIA at the GB, Eg: 
the total energy of a GB system, C: the V/SIA concentration at the GB, £,: the average V-V/SIA-SIA 
binding energy at the GB. 

The GB energy of a pure GB, by definition, is given by: 

Eqs 2 [ EN) - N x E, ]/ Soa» (1) 

The GB energy of the GB with a certain number of defects is given by: 


Eg, -| E WEN Ens E, ]/ S; (2) 


The energy difference is: 


AEG, = Egg, — Egg, = | E (N € n) -Egs(N) TAXE, |! Sop 


-| (N +n) — ES (N) _ 
n 


Fag |S 3) 


-aaen Fam ] 
n 


EEG | C 


The V/SIA formation energy, by definition, is given by: 


By Ea N D-E OT Ea (4) 


EGS(N £ n) - Ec (N) 
In Eq. 3, the term n means the average energy variation induced by the adding of n 


EG (N £1) - E" (N) 


defects. In Eq. 4, the term means the energy variation induced by adding of one defect. 
If there is no interaction among defects, the two terms would be the same. In this case, we obtain: 


_ pGB 
AEG = Esm XC 
5 


(5) 
If there is interaction among defects, then we obtain: 
AEG, ~ (E Sanr - flx E,)xC > (6) 
where flis approximately the fraction of defects pairs relative to the total number of defects. 


From Eqs. 5 and 6, one could see the linear dependence of the GB energy on the defects concentration 


at the GB. 
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